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1.0  INTRODUCTION 


Fluid  flowing  over  a  surface  creates  a  shear  stress  which  acfs  on  that 
surface  according  to  Newton's  Law  of  Viscosity.  If  this  shear  stress,  referred  to 
as  skin  friction,  is  known  over  the  entire  surface,  it  can  be  integrated  over  the 
surface  to  determine  the  skin  friction  component  of  the  total  drag  acting  on  that 
surface.  The  other  major  drag  component,  referred  to  as  pressure  drag  or  form 
drag,  is  due  to  the  boundary  layer  effect  on  the  pressure  distribution  and  also 
depends  on  the  shape  and  orientation  of  a  body  in  a  flow. 

Due  to  the  importance  of  skin  friction  in  determining  the  total  drag  on 
such  surfaces  as  aircraft,  ship  hulls,  etc.,  many  methods  have  been  devised  to 
measure  this  parameter.  The  most  direct  method  is  the  floating  beam  balance 
which  measures  the  force  exerted  by  the  flow  over  a  small  area  of  the  test  section. 
However,  it  is  very  sensitive  to  misalignment  of  the  test  area  with  respect  to  the 
surrounding  surface  and  to  pressure  gradients,  and  is  generally  complicated,  time 
consuming,  and  expensive  to  use.  Other  methods  are  based  on  near-wall  velocity- 
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profile  similarity  assumptions  or  property  analogies.  For  example,  in  the  well 
known  "Law  of  the  Wall",  the  similarity  velocity  profile  scales  on  the  wall  shear 
stress  and  can  be  used  to  determine  the  skin  friction  from  measured  velocity 
profiles.  Similarly,  the  analogy  between  heat  transfer  and  skin  friction  is  often 
exploited  to  provide  a  means  of  measurement  using  heated  metal  films.  These 
are  just  a  few  of  the  more  popular  techniques.  All  of  these  methods,  however, 
are  either  difficult  to  use,  expensive,  intrusive,  or  indirect  and  depend  on 
assumptions  or  analogies.  Thus,  the  search  continued  for  a  direct,  inexpensive, 
nonintrusivc  method  which  can  be  easily  implemented. 

Tanner  and  Blows  (1976)  developed  such  a  method  known  as  Oil  Film 
Laser  Interferometry  based  on  standard  oil  lubrication  theory  and  laser 
interferometry.  Although  the  method  was  both  inexpensive  and  non-intrusive,  it 
was  very  time  consuming.  Monson  and  Higuchi  (1981),  Monson,  Driver,  and 
Szodruch  (1981),  and  Monson  (1984)  simplified  the  method  and  reduced  the  time 
necessary  to  obtain  the  data.  Applying  it  to  a  variety  of  complex  flows  produced 
results  in  good  agreement  with  those  obtained  by  traditional  methods.  They  also 
developed  equations  to  measure  an  unknown  oil  viscosity  by  using  this  method 
with  gravity  as  the  only  force  on  the  oil.  Westphal  (1985,1986)  further  improved 
the  method  by  automating  the  data  reduction  procedure.  In  addition,  he 
investigated  several  parameters  and  set  up  guidelines  for  determining  bad  data 
records.  In  subsequent  studies,  Kim  and  Settles  (1988)  applied  the  method  to 
supersonic  flow  past  a  wedge  and  Cooke  (1988)  applied  the  method  to 
three-dimensional  flow  past  a  wing-body  junction  with  pressure  gradients. 
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This  study  continues  Cooke's  study  by  employing  oil-film  laser 
interferometry  to  investigate  air  flow  around  a  wing-body  junction  and  includes 
an  exploration  of  skin  friction  behind  the  tail  and  in  areas  around  a  line  of  low 
shear.  The  wing  is  composed  of  a  3:2  elliptical  nose  and  a  NACA  0020  tail 
mounted  at  zero  angle  of  attack  (Figure  1).  The  momentum  thickness  Reynolds 
number  for  the  conditions  in  this  study  was  Ree  =  6700  based  on  the  undisturbed 
freestream  velocity  Uref=  21 .1  mis  and  the  momentum  thickness  of  the  approach 
boundary  layer  15.1cm  upstream  of  the  wing  leading  edge.  The  surface  oil  flow 
visualization  in  Figure  2  illustrates  the  primary  features  of  the  flowfield.  The 
flow  includes  a  three-dimensional,  pressure-driven,  turbulent  boundary  layer, 
which  is  symmetric  about  the  wind  tunnel  centerline,  and  is  dominated  by  a 
horseshoe  vortex  formed  around  the  nose  of  the  wing  (Devenport  and  Simpson, 
1988b).  This  vortex  brings  high  momentum  flow  down  along  the  side  of  the  wing, 
into  the  test  wall  boundary  layer  and  out  away  from  the  model.  A  line  of  low 
shear  develops  at  the  outer  edge  of  the  horseshoe  vortex  due  to  the  effect  of  the 
upward  flow  of  the  vortex  at  this  :ocation  (Devenport  and  Simpson,  1988a). 
Figure  3  shows  the  probability  density  functions  of  the  U  component  velocity 
fluctations  measured  in  the  vicinity  of  the  vortex  upstream  of  the  nose  of  the  wing 
at  x'lT  =  —0.20.  The  double-peaked  probability  density  functions  near  the  wall 
indicate  that  the  vortex  is  subject  to  intense,  large-scale,  low-frequency 
unsteadiness.  Exact  details  of  this  unsteadiness  are  currently  unknown.  This 
particular  flow  has  been  the  subject  of  intense  investigation  by  several  researchers 
at  Virginia  Polytechnic  Institute  and  State  University.  Surface  pressure 
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fluctuations,  hot  wire  and  laser  Doppler  anemometer  velocity  measurements,  and 
flow  visualizations  have  been  obtained  (Agarwal  and  Simpson,  1989,  Devenport 
and  Simpson,  1988a,  1988b,  1990). 

In  this  thesis,  the  basic  principles  of  oil  lubrication  theory  and  light 
interference  will  be  discussed  as  they  apply  to  the  measurement  of  skin  friction. 
An  historical  development  of  the  oil  film  laser  interferometer  technique  will 
follow.  The  hardware  and  procedure  used  in  this  study  will  then  be  discussed, 
followed  by  an  analysis  of  the  final  results  and  major  sources  of  error. 
Conclusions  and  recommendations  will  then  be  presented. 
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2.0  THEORY 


Oil-film  laser  interferometry  was  developed  from  two  basic  principles  -  oil 
lubrication  theory  and  laser  interferometry.  Both  principles  were  well  known  and 
accepted  within  the  scientific  community  long  before  Tanner  and  Blows  (1976) 
combined  them  to  measure  skin  friction. 

Oil  lubrication  theory  states  that  a  thin  oil  film  will  flow  over  a  smooth 
surface  under  the  action  of  the  shear  stress  created  by  air  flowing  over  the 
surface.  Squire  (1962)  developed  the  basic  equations  for  oil  flowing  under  the 
action  of  a  shear  stress.  In  a  two  dimensional,  laminar  airflow  with  a  constant 
shear  stress  and  no  shear  gradients,  pressure  gradients,  or  surface  curvature,  the 
oil  film  will  assume  a  linear  profile  (Figure  4).  Starting  with  Newton's  Law  of 
Viscosity 


T 


and  noting  that,  for  the  linear  profile 


(1) 
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du  _  UF  ~  Q 
dy  h  —  0 


(2) 


and 


uF  - 


x-x0 

t 


(3) 


where  t  is  the  time  it  would  take  an  oil  particle  to  travel  from  x0  to  *  at  speed 
uF,  Tanner  and  Blows  (1976)  derived  the  simple  relation  between  the  shear  stress 
and  the  time  dependent  oil-film  thickness  at  a  given  location,  x,  for  such  a  flow 


x-x0 

t  =  u - 

*  h{ts)t 


(4) 


The  measurement  location,  x,  and  the  oil-film  leading  edge,  x0,  are  fixed  for  a 
particular  run  and  can  be  accurately  measured.  The  time  of  flowing,  t,  can  be 
measured  directly  or  calculated  from  the  data  record;  the  latter  method  alleviates 
several  problems  and  will  be  discussed  later.  The  oil  viscosity  as  a  function  of 
temperature  can  be  obtained  from  the  manufacturer's  equation.  Thus,  if  the  time 
dependent  oil-film  thickness  can  be  accurately  determined,  the  time-averaged 
shear  stress  over  the  area  between  x  and  x0  can  be  obtained  from  equation  (4). 

The  time  dependent  oil  film  thickness  can  be  measured  accurately  using  a 
laser  interferometer.  This  is  accomplished  by  directing  a  laser  beam  onto  an  oil 
film  flowing  over  a  partially  reflective  glass  surface  as  depicted  in  Figure  5.  Part 
of  the  laser  beam  is  specularly  reflected  from  the  air-oil  interface  B,  resulting  in 
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beam  BE,  and  the  remainder  of  the  beam  is  transmitted  into  the  oil  film  as  beam 
BC.  A  similar  reaction  occurs  at  the  oil-surface  interface  C,  creating  beams  CD 
and  CJ.  At  the  oil-air  interface,  part  of  beam  CD  is  transmitted  as  beam  DF  and 
part  is  reflected.  This  process  will  continue,  but  due  to  losses  from  the  multiple 
reflections,  beams  other  than  BE  and  DF  can  be  neglected.  Since  BE  and  DF 
are  parallel  and  coherent,  they  interfere  with  each  other,  and  the  resulting 
interference  pattern  is  detected  by  the  receiving  optics.  As  the  oil  film  thins  due 
to  the  action  of  the  shear  stress  acting  upon  it,  the  differences  in  path  length 
between  beams  BE  and  DF  decrease.  This  path  length  difference  changes  the 
relative  phases  of  beams  BE  and  DF  at  parallel  wavefronts  causing  the  intensity 
of  the  combined  beams  to  vary  as  a  sinusoid  with  a  decaying  period  (Figure  6). 
As  the  oil  film  continues  to  thin,  the  path  length  difference  goes  through  full 
interference  cycles  -  varying  from  maximum  constructive  interference  through 
maximum  destructive  interference  and  back  to  maximum  constructive 
interference.  This  cycle,  referred  to  as  one  "fringe,"  corresponds  to  a  change  in 
the  path  length  difference  of  one  laser  wavelength.  Note  that,  due  to  a  difference 
in  indices  of  refraction  between  the  air  and  the  oil,  the  wavelength  of  the  laser  in 
the  oil  differs  from  the  wavelength  in  air  according  to  the  relation 


When  light  travelling  in  a  medium  with  refractive  index  /?,  encounters  an 
interface  with  a  material  of  refractive  index  n2  where  n2  <  n,  both  transmission 
and  reflection  occur,  but  the  phase  of  the  reflected  beam  is  shifted  180°.  If 


THEORY 


7 


n2  >  n„  then  no  phase  shift  occurs.  Note  that  the  phase  of  the  transmitted  beam 
is  never  shifted,  regardless  of  nx  and  n2. 

None  of  the  reflections  which  result  in  beams  BE  and  DF  cause  a  phase 
change.  Therefore,  when  the  difference  in  path  length  between  beams  BE  and 
DF  is  an  integral  number  of  wavelengths,  a  constructive  maximum  will  occur. 
The  path  length  difference,  A/,  is  (refer  to  Figure  5  for  geometry  and  notation) 

A/  =  n0(BCD)  -  BI  (6) 

However,  ID  and  BH  represent  successive  positions  of  the  same  wavefront 
(Jenkins  and  White,  1957).  Therefore,  the  beams  HD  and  Bl  contain  the  same 
number  of  wavelengths 


n0(HD)  =  BI  (7) 

and  the  path  length  difference  becomes 

A/  =  n0(BCH)  +  n0{HD)  -  BI  =  n(BCH)  (8) 

By  extending  DC  until  it  intersects  a  vertical  line  through  B,  and  noting  that 
GC  =  BC  due  to  equal  angles  of  incidence  and  reflection  at  point  C,  the  Final  form 
of  the  path  length  is 

A/  =  n0(GC  +  CH)  =  n0{GH)  =  n0\2h  cos(<?0)]  (9) 

As  stated  previously,  a  maxima  will  occur  when  the  path  length  difference  is  an 
integral  multiple  of  wavelengths 
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mX  =  2n0hcosd0  m—  1,2,3,...  (10) 

Thus,  a  change  in  the  path  length  of  A0  corresponds  to  a  change  in  oil  film 
thickness  of 


2 n0  cos  60 


(11) 


where 


e0 


sin  0( 

no  " 


(12) 


This  change  in  height  is  referred  to  as  a  thickness  change  of  one  "fringe." 
By  recording  the  interference  level  as  a  function  of  time  and  allowing  the  oil  film 
to  thin  to  less  than  one  fringe  in  height,  the  height  of  the  oil  film  at  any  time  can 
be  computed  by  simply  counting  the  fringes  backward  (peaks  labeled  as  "N  =" 
in  Figure  6)  beginning  with  the  final  fringe,  to  the  desired  time.  The  height  of  the 
oil  film  at  that  time  is  simply  the  number  of  interference  fringes,  N  (referred  to 
as  the  fringe  number),  times  the  "fringe  thickness,"  hx  (h (r„  x)  =  hxNt).  Thus,  the 
oil  lubrication  equation  becomes 


T  = 


x-x0 
*  hxNt 


(13) 


Since,  hx,  x  -  n,  and  r  are  constant  for  a  given  test  run,  the  product  of 
fringe  number,  N,  and  time,  r,  must  also  be  a  constant,  C{x),  at  a  given  location. 
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C(  x)  =  constant  —  Nt 


(14) 


Note  that  most  interferometers  reflect  one  of  the  interfering  beams  off  of 
a  stationary  reference  and  measure  the  change  in  distance  of  the  test  surface  often 
to  less  than  one  micron.  Thus,  interferometers  are  generally  very  sensitive  to 
vibration  of  the  test  surface.  However,  the  system  used  to  measure  thinning  of 
an  oil  film  is  relatively  insensitive  to  vibrations  of  the  test  surface  since  the 
interference  occurs  due  to  the  distance  between  the  oil-air  interface  and  the 
oil-surface  interface.  Thus,  if  the  test  surface  is  shifted  upward  a  small  amount 
during  vibration,  the  oil  film  will  also  be  shifted  upward  by  the  same  amount. 
While  this  changes  the  actual  distance  travelled  by  the  laser  beams  from  the  laser 
to  the  receiving  optics,  it  will  not  affect  the  path  length  difference  which  is  what 
causes  the  interference  to  occur. 
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3.0  HISTORICAL  DEVELOPMENT 


The  first  oil  film  laser  interferometer,  designed  and  built  by  Tanner  and 
Blows  (1976),  consisted  of  a  He-Ne  laser  beam  directed  through  a  beamsplitter 
to  create  two  laser  beams  a  fixed  distance  apart;  this  distance  was  measured 
directly.  A  drop  of  oil  was  then  placed  on  the  test  section  and  the  beams  were 
directed  onto  the  test  section  where  the  upstream  beam  was  visually  positioned 
at  the  oil  drop's  leading  edge,  x0.  Thus,  only  the  downstream  beam  was  used  to 
measure  the  thinning  of  the  oil  film.  For  this  reason,  it  is  referred  to  as  a  single 
beam  system.  The  photodetector  output  from  the  receiving  optics  was  initially 
monitored  on  a  strip  chart  recorder,  and  the  positively  sloped  zero  crossings  were 
read  off  manually  to  denote  the  interference  fringes.  Note  that,  using  the 
positively  sloped  zero  crossings  to  identify  the  fringes  yields  fringe  numbers  which 
will  be  integers  plus  1/4  (i.e.  1.25,  2.25,  etc.)  as  a  result  of  equation  (10).  Later, 
an  analog  fringe-detection  circuit  and  computer  were  used  to  identify  the  zero 
crossings. 
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Two  separate  data  reduction  methods  were  applied  to  this  single  beam 
system.  First,  the  oil  was  allowed  to  thin  to  less  than  one  fringe,  as  mentioned 
earlier,  and  an  effective  time  of  flowing  was  computed  using  a  best  fit  to  equation 
(14).  (For  clarity,  when  the  time  or  fringe  number  is  calculated  instead  of  being 
measured  directly,  an  "eff"  subscript  will  be  used  to  denote  an  effective  quantity 
as  opposed  to  a  measured  quantity.)  This  procedure  resulted  in  very  long  test 
runs-some  as  long  as  60  minutes.  As  an  alternative,  the  time  of  flowing  was 
referenced  to  the  time  the  wind  tunnel  was  started,  and  the  effective  fringe 
numbers,  we re  calculated  from  equation  (14).  This  allowed  for  shorter  test 

runs  since  the  oil  did  not  have  to  thin  to  less  than  one  fringe.  However,  this 
method  underpredicts  the  actual  time  of  flowing  due  to  the  initial  slope  of  the  oil 
film.  This  time  discrepancy  was  shown  to  be  where  a  is  the  initial  slope  of 
the  oil  film  leading  edge  (Tanner  and  Blows,  1981),  and  was  computed  to  be 
typically  between  2  seconds  and  5  seconds  for  the  specific  flow  conditions  studied. 

Tanner's  method  had  several  drawbacks  that  prevented  its  widespread  use. 
First,  it  often  entailed  very  long  test  runs  —  especially  for  low  shear  flows. 
Second,  visually  locating  the  oil  film  leading  edge  was  found  to  be  a  difficult  task 
in  practice  and  the  leading  edge  was  also  subject  to  move  slightly  due  to  oil 
spreading  after  it  had  been  located  with  the  upstream  laser  beam  and  before  the 
wind  tunnel  was  started.  Furthermore,  the  effective  flow  time  of  the  oil  differs 
from  the  length  of  time  the  wind  tunnel  is  operating  due  to  transient  effects  and 
the  initial  slope  of  the  leading  edge  of  the  oil  drop. 
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Monson  and  Higuchi  (1981)  were  able  to  eliminate  the  major  problems 
with  Tanner's  system  by  positioning  both  beams  downstream  of  the  leading  edge 
and  using  both  to  record  data.  Monson  perpendicularly  polarized  the  two  beams 
and  separated  the  signals  they  carried  with  polarizing  filters.  Since  both  beams 
were  used  to  obtain  data,  this  is  referred  to  as  a  dual-beam  system.  The  time  of 
the  signal  peaks,  instead  of  the  positively  sloped  zero  crossings  used  by  Tanner 
and  Blows  (1976),  were  used  to  denote  the  fringe  numbers  to  eliminate  difficulties 
arising  from  any  drift  in  the  mean  signal  level.  Therefore,  Monson  and  Higuchi's 
fringe  numbers  were  integers.  Noting  that  t  (taken  as  an  average  over  the  area 
between  the  two  beams),  /i,  and  x0  are  identical  for  both  beams,  the  leading  edge 
position  can  be  eliminated  algebraically  by  solving  one  of  the  following  equations 
for  .x0  and  substituting  it  into  the  other. 


t  =  m 


*i  -*o 


(15) 


T  -ll 


*2~*0 

hflefptefp 


(16) 


After  simplification,  Monson's  equation  for  a  dual-beam  laser 
interferometer  becomes 


*2~*1 

hANef alefJ2  ~  Neff\teff[) 


(17) 
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Thus,  the  distance  between  the  two  laser  beams  is  required,  but  the  leading 
edge  need  not  be  located.  Also,  Monson  ran  each  individual  test  run  only  long 
enough  to  identify  20  "good"  fringes  (smooth,  easily  identifiable,  etc.)  on  the 
downstream  beam,  recording  his  output  on  a  strip-chart  recorder.  He  then 
developed  equations  to  calculate  the  effective  time  based  upon  the  elapsed  time 
for  20  fringes  on  the  downstream  beam  to  occur,  and  for  the  fringe  numbers  for 
each  beam  based  on  the  elapsed  time  for  10  fringes  from  each  beam  to  occur 
(Monson  and  Higuchi,  1981).  This  considerably  reduced  the  time  necessary  to 
obtain  the  data,  eliminated  all  the  problems  related  to  locating  the  leading  edge 
of  the  oil  film,  and  allowed  the  time  origin  to  be  computed  without  having  to  let 
the  oil  film  thin  to  less  than  one  fringe. 

Monson  and  coworkers  applied  their  dual  beam  system  to  a  wide  variety 
of  flows  and  derived  equations  and  correction  terms  for  various  flow  conditions. 
They  showed  that  if  the  direction  of  measurement  (denoted  by  an  imaginary  line 
connecting  the  two  laser  beams)  were  different  than  the  actual  direction  of  the 
skin  friction,  the  interferometer  system  would  measure  the  component  of  skin 
friction  in  the  direction  of  measurement.  In  addition,  they  derived  first  order 
correction  terms  for  shear  gradients,  pressure  gradients,  linear  temperature 
gradients,  and  body  forces.  Up  to  this  point  in  the  development  of  the  Oil  Film 
Laser  Interferometer,  'bad'  data  records  were  rejected  subjectively  upon  visual 
inspection  of  the  strip-chart  recorder  output.  Uneven  signal  visibility  and 
obvious  noisiness  were  considered  the  main  reason  for  rejection  of  bad  data 
records. 
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Murphy  and  Westphal  (1985)  and  Westphal  and  Bachalo  (1986)  were  the 
next  to  improve  the  system.  They  obtained  a  more  equal  signal  visibility  between 
the  two  beams  than  Monson  by  spatially  resolving  the  laser  beams.  However, 
Westphal's  most  important  contribution  was  to  automate  the  data  reduction 
procedure.  Instead  of  using  a  strip  chart  recorder,  he  converted  the  photodetector 
output,  using  an  analog-to-digital  converter,  to  a  digital  signal  which  was  stored 
directly  on  a  computer  hard  disk.  He  wrote  FORTRAN  routines  to  identify 
fringes  from  the  data  record  and  a  tw’o-variable  least  squares  error  minimization 
routine  to  obtain  the  best  fit  to  equation  (14).  He  defined  the  effective  time  and 
fringe  numbers  as 


II 

1 

(18) 

V/  ~  r0  +  Oc 

(19) 

where  k  is  the  chronological  fringe  number  within  the  data  record  (see  Figure  6), 
f*  is  the  time  of  the  k'h  positively  sloped  zero  crossing  referenced  to  the  beginning 
of  the  data  record,  N0  is  a  constant  to  give  the  actual  height,  in  hx  s,  and  t0  is  the 
time  origin  or  effective  time  of  oil  flowing.  Note  that,  if  transient  effects  and  the 
initial  slope  of  the  oil  film  leading  edge  are  negligible,  t0  is  identical  to  the  time  the 
wind  tunnel  was  started.  The  error  minimization  routine  iterates  on  N0  and  f0  to 
find  the  minimum  normalized  rms  error  of  the  constant  C{x) 

7)  -  (c)2 
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It  is  important  to  note  that  previously,  Ntfr  and  ttfJ  were  calculated  using 
only  two  or  three  distinct,  widely  spaced  fringes.  Using  the  power  and  speed  of 
the  computer,  Westphal  used  all  of  the  fringes  to  calculate  NtfJ  and  te}r  In 
addition,  Westphal  computed  a  time  origin  t0  from  each  data  record,  thus 
producing  two  values  for  the  time  origin.  Since  the  time  origin  is  a  property  of 
the  entire  oil  flow,  it  is  independent  of  the  measurement  location  fx,  or  x2). 
Realizing  this,  Westphal  determined  that  if  the  difference  exceeded  a  certain  level, 
then  at  least  one  of  the  data  records  was  "bad."  The  criterion  he  established  for 
rejection  was 


I  frl  ~  (Q2  1 
v/ at  k  -  1) 


>  .02 


(21) 


Westphal  also  established  a  maximum  allowable  error  on  the  constant 
C{x )  for  each  beam,  defined  in  equation  (20).  This  criterion  limited  the  error  to 
e  <  .004  ;  if  either  beam  failed  to  meet  this  requirement,  the  data  record  was 
considered  unacceptable.  In  practice,  he  found  most  good  runs  to  have  e<.003. 
These  rejection  criteria  were  based  upon  a  maximum  allowable  error  of  ±  5%  on 
the  final  skin  friction  result  for  a  single  run.  This  accounted  for  variations  in  oil 
viscosity,  beam  separation  measurements  (x2  -  uc,),  etc.  These  criteria  will 
produce  a  result  within  a  ±  5%  uncertainty  only  for  a  two  dimensional,  laminar 
flow  with  no  pressure  or  shear  gradients  and  no  gravitational  effects  such  as  the 
flow  Westphal  studied. 

It  is  important  to  note  here  that  Westphal  (1986)  defines  the  variable  X0 
to  be  one  wavelength  in  the  oil  (which  is  correctly  defined  in  equation  (5)),  but 
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then  gives  an  equation  for  a  thickness  variation  (referred  to  here  as  hx)  of  the  oil 
film  corresponding  to  a  change  of  one  X0  and  uses  the  same  variable  name,  X0,  for 
this  change  in  height.  However,  his  subsequent  analysis  is  consistent  with  the 
definition  of  Xe  as  the  change  in  height.  For  clarity,  hx  is  used  in  this  study. 

Cooke  (1988)  applied  Oil- Film  Laser  Interferometry  to  flow  past  a  wing 
body  junction.  He  used  the  programs  written  by  Westphal  (1986)  combined  with 
a  linear  pressure  gradient  correction  term  developed  by  Monson,  Driver,  and 
Szodruch  (1981).  He  determined  the  skin  friction  as  if  no  pressure  gradient 
existed,  and  then  modified  the  result  as  follows 


t'  = 


(22) 


where 


£  = 


WeffB 
2n0r  cos  60 


dp_ 

dx 


(23) 


Murphy  and  Westphal  (1985)  showed  that  for  sufficiently  thin  oil  films  (less  than 
20  microns),  this  system  would  measure  the  time-averaged  shear  stress  in  a 
three-dimensional  flow.  For  oil  films  thicker  than  20  microns,  they  showed  that 
surface  waves  may  appear  at  the  oil-air  interface  depending  on  the  particular 
flow  conditions. 

Cooke  did  not  actually  measure  both  the  magnitude  and  direction  of  this 
shear  stress  in  his  experiments.  Instead,  he  pre-determined  the  direction  of  flow 
from  surface  oil  flow  patterns  produced  by  earlier  researchers  working  on  the 
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same  flow.  Knowing  the  flow  angle,  he  aligned  the  laser  beams  with  the  direction 
of  flow  and  measured  the  skin  friction  magnitude.  He  also  experimented  with 
various  scanning  devices  which  could  take  data  at  several  closely  spaced  data 
points,  but  did  not  find  an  acceptable  device.  He  also  experimented  with 
measuring  the  oil-film  thickness  from  below  the  test  section.  This  requires  the 
thickness  of  the  glass  plate  to  be  measured  and  subtracted  from  the  thickness 
measured  with  the  oil  on  the  test  section.  His  final  decision  was  to  use  a 
dual-beam  system  mounted  above  the  test  section.  The  system  was  placed  on  a 
wooden  frame  which  surrounded  the  wind  tunnel  but  was  not  in  contact  with  it. 
This  eliminated  the  possibility  of  transmitting  vibrations  from  the  wind  tunnel 
motor  to  the  interferometer  system. 

With  the  exception  of  Cooke,  this  method  had  only  been  applied  to  well 
understood  flows  such  as  two-dimensional  flows,  axisymmetric  flows,  flows  with 
one-dimensional  pressure  gradients,  etc.  Cooke  applied  this  method  to  a  complex 
three-dimensional  flow,  but  relied  upon  surface  oil  flow  visualizations  to 
determine  the  flow  angles  and  then  aligned  the  laser  beams  with  the  flow. 

For  the  situation  in  which  the  flow  angle  is  unknown,  Monson,  Driver, 
and  Szodruch  (1981)  derived  an  important  result  for  two  dimensional  flow. 
Figure  7  illustrates  the  geometry  and  notation  used  to  derive  the  flow  angle  and 
shear  stress  for  this  situation.  For  a  constant  shear  stress  and  constant  oil 
viscosity  (i.e.  constant  temperature),  the  oil  lubrication  equations  at  points 
(x'„  z\)  and  (jc'2,  z'2)  (refer  to  Figure  7  for  notation  and  geometry)  are 
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Ms\  ~  ^oi) 

T“  Vi 


(24) 


His2  ~  502) 

T  “  h2t2 


(25) 


Since  the  flow  direction  is  unknown,  the  streamwise  variables,  s0 ,  and  s2  must 
be  solved  for  in  terms  of  the  measurement  direction  variables  x0,  x,,  and  x2 


x\  ~  xo 

sl  ~  so\  =  cos  y  ^ 

x2  ~  x0 

s2  -  s02  =  cosy  {-11> 

Substituting  into  eqautions  (21)  and  (22)  and  eliminating  x0  as  previously  done 
yields 


t  = 


h2t2-hxtx  cosy 


(28) 


Thus,  the  total  shear  t  is  unknown  since  y  is  unknown,  but  t,  =  x  cos  y  which 
gives  the  final  result 


_  m(x2~xi) 

x  h2t2-h\tx 


(29) 


Thus,  oil-film  laser  interferometry  measures  the  component  of  the  shear  stress  in 
the  direction  of  an  imaginary  line  connecting  the  laser  beams  on  the  test  surface 
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(known  as  the  measurement  direction).  It  is  very  important  that  the  initial  line 
of  oil  be  placed  perpendicular  to  the  measurement  direction  as  will  be  shown  in 
Chapter  7. 

In  this  thesis,  measurements  were  initially  taken  in  the  x  and  -z 
directions  (using  the  nomenclature  of  Figure  1)  so  the  measurements  in  each 
direction  could  be  combined  vectorially  to  obtain  the  magnitude  and  direction  of 
the  shear  stress  at  each  point.  However,  when  the  angle  between  the  flow  and  the 
x-axis  was  less  than  approximately  20°  and  the  laser  interferometer  was  aligned 
in  the  -z  direction,  the  oil  did  not  flow  over  both  data  beams  (see  Figure  8). 
Therefore,  data  could  not  be  obtained  in  the  —  z  direction,  and  the  angles  and 
magnitudes  of  the  actual  shear  stress  could  not  be  determined. 

Since  both  principal  axis  directions  could  not  be  used,  measurements  were 
taken  in  three  directions  far  enough  apart  to  provide  acceptable  spatial 
resolution.  Data  were  obtained  at  0°,  —45°,  and  -60°  with  respect  to  the  wind 
tunnel  axis.  Note  that  measurements  in  only  two  directions  are  sufficient  to 
obtain  the  skin  friction  magnitude  and  direction.  Measurements  in  the  third 
direction  were  included  to  provide  redundancy.  The  rejection  criteria  developed 
by  Westphal  ( e  <  .004)  were  applied  to  this  data  and  most  good  runs  were 
observed  to  have  e  <  .003  for  each  beam,  in  agreement  with  his  findings. 
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4.0  APPARATUS 


All  the  experiments  from  this  study  were  done  in  the  VP1&SU  Low-Speed 
Boundary-Layer  Wind  Tunnel  (Figure  9).  The  tunnel  is  an  open  circuit  tunnel 
with  a  rectangular  test  section  8m  long  x  .91m  wide  x  .26m  high  (over  the  central 
portion).  Air  is  taken  into  a  centrifugal  blower  through  a  Filter.  Upstream  of  the 
test  section  is  a  fixed-setting  damper,  a  plenum,  a  honeycomb  section,  seven 
screens,  and  a  4:1  contraction  ratio  nozzle  which  leads  to  the  test  section.  In  the 
absence  of  the  wing,  a  nearly  zero  pressure  gradient  exists  in  the  test  section.  At 
the  inlet  to  the  test  section,  the  flow  is  uniform  to  within  .5%  in  the  spanwise 
direction  and  1%  in  the  vertical  direction  with  a  .2%  turbulence  intensity  at 
27 m/s.  The  flow  is  tripped  as  it  enters  the  test  section  by  a  .63cm  step,  and  is 
accelerated  to  test  speed  by  a  further  1.5:1  contraction  produced  by  the  shape  of 
the  upper  wall.  The  wind  tunnel  is  the  same  as  that  used  by  Cooke  (1988). 

The  upper  wall  of  the  wind  tunnel  is  made  of  plexiglas  reinforced  with 
aluminum  channel.  The  side  tunnel  walls  are  made  of  glass  and  lined  internally 
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with  6.4mm  thick  plexiglas  plates.  When  the  wing  is  mounted  in  the  test  section, 
the  plexiglas  liners  near  the  wing  are  removed  to  minimize  the  effects  of  blockage 
induced  pressure  gradients.  This  effectively  widens  the  test  section  by  12.7mm  in 
the  area  between  a  location  330mm  upstream  of  the  leading  edge  of  the  wing  and 
203mm  downstream  of  the  trailing  edge  of  the  wing.  These  distances  were  chosen 
to  make  the  side  walls  follow,  very  approximately,  a  streamline  produced  by  the 
wing  in  an  unbounded  flow.  The  corners  produced  by  the  removal  of  the 
plexiglas  liners  are  taped  over  adhesive  tape  to  reduce  added  turbulence  and 
recirculation  due  to  the  surface  discontinuity.  The  streamlines  followed  by  the 
side  walls  were  calculated  using  two-dimensional  potential  flow  theory. 

The  floor  of  the  test  section  is  made  of  19mm  thick  fin-form  plywood.  The 
area  immediately  surrounding  the  permanent  wing  mount  has  been  cut  away  so 
different  materials  may  by  used  for  the  actual  test  surface.  The  test  area  inserts 
used  in  this  study  consist  of  two  parts.  The  first  is  made  of  3/4"  Finished  plywood 
secured  flush  with  the  surrounding  floor  by  countersunk  bolts  located  at  the 
edges.  A  cutout  section  with  a  1/2"  lip  holds  the  second  piece  -  a  1/4"  thick 
Denton  Vacuum  MLBS-30  partially  reflective  glass  plate.  The  glass  plate,  chosen 
by  Cooke  to  match  the  intensity  of  the  interfering  beams,  has  a  30%  reflective 
coating  on  the  test  side.  Cooke  took  data  only  in  front  of  the  wing,  using  the 
plywood  insert  shown  in  Figure  10(a)  The  present  study  investigated  positions  in 
front,  along  the  side,  and  behind  the  wing.  Therefore,  two  additional  plywood 
inserts  were  constructed,  with  cutouts  as  shown  in  Figures  10(b)  an  ;  10(c). 
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The  wing  used  in  this  study  consists  of  a  3:2  elliptical  nose  upstream  of  the 
maximum  thickness  and  a  NACA  0020  tail.  The  maximum  thickness  is  7.17  cm, 
the  chord  is  30.48  cm,  and  the  height  of  the  wing  is  22.9  cm.  This  height  allowed 
for  a  3.7  cm  gap  between  the  top  of  the  wing  and  the  wind  tunnel  roof  to  prevent 
the  formation  of  a  junction  vortex  at  the  top  of  the  wing.  Such  a  vortex  could 
have  interfered  with  the  flow  further  downstream  on  the  test  wall.  Sandpaper 
strips  (120  grade)  located  10  mm  upstream  of  the  maximum  thickness  ensured  a 
constant  turbulence  transition  point  on  the  wing  surface. 

The  plywood  insert  and  glass  plate  were  secured  flush  with  the 
surrounding  tunnel  floor  and  the  wing  was  bolted  to  a  permanent  mount  in  the 
tunnel  floor.  The  wing  was  mounted  with  its  leading  edge  1 .39m  downstream  of 
the  tunnel  throat  and  its  chord  parallel  to  the  tunnel  centerline.  Small  gaps 
between  the  inserts  and  the  surrounding  floor  were  covered  with  Scotch  brand 
"Magic  Mending"  tape  to  provide  a  smooth  surface  continuity. 

A  wooden  frame  was  constructed  around  the  wind  tunnel  to  support  the 
laser  assembly.  On  top  of  this  structure,  Cooke  secured  the  laser  assembly  such 
that  it  was  directly  above  the  test  location.  In  order  to  provide  access  to  the 
additional  downstream  locations  investigated  in  this  study,  Cooke's  structure  was 
extended  two  and  a  half  feet  in  the  streamwise  direction  (see  Figure  1 1  for  the 
final  structure).  With  Cooke's  setup,  the  time  required  to  move  between  locations 
was  about  45  minutes.  This  large  relocation  time  occurred  because  the  laser 
assembly  had  to  be  physically  picked  up  and  moved,  carefully  re-aligned  with  the 
flow  direction,  and  clamped  in  place.  To  reduce  this  relocation  time,  a  movable 
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platform  was  mounted  on  top  of  the  wooden  structure  on  a  rail  system  (Figure 
12).  The  movable  platform  could  slide  along  the  rails  in  the  streamwise  (top  set 
of  rails)  and  spanwise  (bottom  set  of  rails)  directions,  thus  allowing  the  laser 
assembly  mounted  on  top  of  the  movable  platform  to  be  relocated  without 
changing  the  alignment  angle.  This  was  done  by  attaching  two  parallel  metal 
angles  to  the  bottom  of  the  platform.  These  angles  subsequently  fit  into  an 
additional  set  of  oppositely  oriented  angles  with  the  same  spacing,  forming  a  set 
of  rails.  Another,  orthogonally  oriented  pair  of  angles  was  bolted  to  the  bottom 
of  the  first  set  of  rails  and  fit  into  a  pair  of  angles  rigidly  mounted  atop  the 
wooden  frame. 

To  allow  for  easy  rotation  of  the  laser  assembly  with  respect  to  the 
platform,  a  1/2"  thick  "U-shaped"  aluminum  piece  was  attached  to  the  laser 
assembly  base  plate.  One  inch  holes  were  drilled  in  both  the  platform  and  the 
U-shaped  fixture,  and  a  1"  bolt  through  both  holes  firmly  attached  the  laser 
assembly  to  the  movable  platform.  The  laser  beams  were  aligned  with  the  three 
different  measurement  directions  in  a  similar  manner.  Three  1/4"  holes  in  the 
platform,  when  aligned  with  a  1/4"  hole  in  the  laser  assembly  base  plate,  oriented 
the  laser  beams  at  0°  ,  45°,  or  60°  with  respect  to  the  wind  tunnel  centerline. 
Wheels  mounted  to  the  underside  of  the  laser  assembly  base  plate  further 
simplified  the  rotation. 

The  platform  was  made  of  3/4"  birch  plywood  with  a  1/4"  aluminum  plate 
affixed  to  the  corner  with  a  1"  bolt  to  provide  additional  strength.  Two  90°  arcs 
cut  through  the  platform  and  aluminum  plate  gave  the  laser  beams  clear  access 
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to  the  test  section  at  any  angle  between  0°  and  90°  with  respect  to  the  tunnel 
centerline.  Figure  13  shows  the  movable  platform  mounted  on  the  rail  system 
above  the  wooden  frame.  Figure  14  shows  the  author  applying  oil  to  the  test 
surface  and  depicts  the  relative  positions  of  the  wooden  frame,  rail  system, 
movable  platform,  and  laser  assembly  with  respect  to  the  wind  tunnel. 

The  basic  laser  interferometer  system  used  was  that  developed  by  Cooke. 
Figure  15  shows  the  laser  assembly  from  the  side  and  Figure  16  shows  the  laser 
assembly  looking  down  the  optical  rail  axis.  Two  parallel  optical  rails  were  firmly 
mounted  on  a  common  1/4"  aluminum  base  plate,  and  all  optical  components 
were  mounted  and  aligned  on  these  rails.  A  Spectra  Physics  Model  105-1  5  mW 
He-Ne  laser  (A  =  6328A)  was  mounted  such  that  lateral  and  vertical  adjustments 
could  easily  be  made.  The  laser  beam  paths  shown  in  Figure  17  are  described 
as  follows.  The  laser  beam  first  passes  through  a  beam  expander  and  focusing 
assembly  consisting  of  a  20X  objective  lens  and  an  output  lens  with  a  focal  length 
f  =  301mm.  This  assembly  is  mounted  on  the  same  optical  rail  as  the  laser  beam. 
The  beam  is  then  reflected  from  a  plane  mirror,  mounted  on  the  same  optical  rail 
as  the  laser,  to  the  beamsplitter,  a  6mm  thick  optical  window,  which  is  mounted 
on  the  second  optical  rail.  The  two  beams  are  then  directed  down  to  the  test 
surface  via  an  adjustable  plane  mirror.  Although  the  beamsplitter  creates  a 
pathlength  difference  which  prevents  both  beams  from  being  simultaneously 
focused  on  the  test  surface,  this  difference  is  not  crucial.  Rather,  focal  lengths 
were  adjusted  such  that  the  laser  beams  have  approximately  the  same  diameter 
at  the  test  surface  -  one  focused  slightly  above  the  test  surface  and  the  other 


APPARATUS 


25 


focused  the  same  distance  below  it.  The  returning  beams  are  reflected  from  a 
plane  mirror  similar  to  the  one  used  for  locating  the  beams  on  the  test  surface, 
and  are  then  re-focused  through  a  convex  lens  with  focal  length  f  =  401mm. 
These  focused  beams  are  split  spatially  by  two  plane  mirrors  mounted  90  degrees 
apart,  and  are  again  reflectd  to  finally  impinge  on  United  Detector  Technology 
silicon  PiN-lODP  photodiodes  operating  in  photovoltaic  mode.  The  photodiode 
amplifier  circuits  were  the  same  as  those  used  by  Cooke  (1988).  The  spot  size  of 
the  laser  beams  on  the  photodiodes  was  controlled  by  an  adjustable  iris.  The  iris 
opening  used  in  these  experiments  was  approximately  1.5mm  to  2mm. 

Each  photodiode  signal  was  amplified  to  approximately  2.5  volts  and 
subsequently  offset  to  allow  bi-polar  analog-to-digital  conversion.  This  offset 
adjustment  was  calibrated  through  use  of  a  voPmeter  while  in  complete  darkness 
(lights  and  laser  turned  off).  The  voltmeter  responded  too  slowly  to  be  useful 
during  data  acquisition.  A  Krohn-Hite  model  3202  fourth  order  Simple  RC  filter 
was  operated  in  the  low-pass  mode  at  a  cutoff  frequency  of  20  Hz  to  eliminate 
noise  from  vibration  of  the  tunnel  motor,  etc.  The  signal  was  monitored  in  real 
time  on  an  oscilloscope;  data  acquisition  was  initiated  when  interference 
variations  were  detected.  A  Data  Translation  DT  2801 -A  analog-to-digital 
converter  board  installed  in  an  IBM  PC  XT  operating  at  10  Hz  recorded  the 
digital  signal  to  a  computer  floppy  disk  for  later  analysis.  Figure  18  shows  a 
block  diagram  of  the  components  as  they  were  set  up  for  this  study. 

The  oil  used  was  Dow  Corning  200  fluid  (manufactured  by  Dow  Corning 
Corporation,  Midland.  MI),  a  50  centi-Stoke  nominal  viscosity  silicone  oil  with  a 
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nominal  temperature  of  25°  Celsius.  The  oil  was  chosen  based  on  Cooke's 
recommendation  for  the  shear  levels  involved.  The  temperature  was  monitored 
from  a  thermometer  attached  to  the  inside  of  the  test  section  approximately  two 
feet  downstream  of  the  tail  of  the  model.  The  manufacturer's  equation  for  oil 
viscosity-temperature  dependence  was  used  to  determine  the  oil  viscosity  for  each 
run 


KT0)  = 


_ ^ nom _ 

e.0Ub(To-Tnom) 


(30) 
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5.0  PROCEDURE 


A  transparency  with  the  measurement  locations  was  attached  to  the 
bottom  of  the  glass  plate  to  assist  in  positioning  the  laser  beams  in  the  proper 
location  (Figure  19).  The  laser  beams  were  positioned  such  that  the  measurement 
location  (marked  on  the  transparency)  was  directly  between  the  two,  since  the 
average  shear  stress  between  the  two  beams  is  actually  measured.  The  laser  and 
the  lights  were  turned  off  and  the  zero  level  for  the  photodetectors  was  set  to 
-1.25  volts.  The  laser  was  then  turned  on  and  fine  adjustments  were  made  to  the 
receiving  optics  alignment. 

The  beam  spacing  was  then  measured  directly  by  placing  a  piece  of  paper 
on  the  test  section  and  measuring  the  distance  between  the  upstream  edges  of  the 
laser  beams  with  a  caliper.  These  measurements  were  accurate  to  .05  mm.  The 
beam  spacing  used  in  these  experiments  was  4.35mm.  The  test  section  and 
surrounding  area  were  then  cleaned  using  commercial  glass  cleaner  and  lens 
paper.  The  oil  was  kept  in  a  syringe  to  protect  it  from  dust  while  still  allowing 
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easy  application.  The  oil  was  applied  from  the  syringe  to  the  edge  of  a  taped 
razor  blade  (the  tape  was  to  prevent  scratching  the  reflective  surface  of  the  glass 
plate).  The  razor  blade  was  then  pressed  down  on  the  test  section  between  one 
and  two  beam  spacings  ahead  of  the  upstream  beam,  creating  a  line  of  oil 
perpendicular  to  the  measurement  direction.  The  tunnel  was  then  started.  Once 
the  tunnel  reached  operating  conditions,  the  dynamic  pressure  (qx  =  ^  pU}tj)  was 
measured  using  a  pitot-static  tube  mounted  in  the  wind  tunnel  throat.  The 
temperature  was  also  recorded;  it  was  normally  between  2 4°C  and  26°C.  When 
the  intensity  of  the  downstream  beam,  viewed  in  real  time  on  an  oscilloscope, 
reached  a  smooth  oscillation,  the  lights  in  the  lab  were  turned  off  (to  prevent  the 
photodetectors  from  picking  up  stray  light)  and  the  computer  program, 
OILAD32  (see  Appendix  A)  ,  which  recorded  the  data  to  a  computer  disk  was 
begun.  The  amount  of  time  between  when  the  tunnel  reached  the  designated 
dynamic  pressure  and  when  the  beginning  of  the  data  recording  was  typically 
between  one  and  two  minutes.  When  the  data  was  finished  recording,  the  lights 
were  turned  on  and  the  initial  temperature  was  verified. 

Figures  20,  21,  22,  and  23  are  examples  of  typical  raw  data  records 
obtained  using  this  method.  Figures  24  and  25  (upstream  and  downstream  data 
records  for  the  same  test  run)  show  an  example  of  a  data  record  that  cannot  be 
reduced.  The  most  likely  reason  for  the  signal  degradation  in  Figure  25  is  that 
the  receiving  optics  gradually  shifted  over  the  test  run  due  to  vibration. 
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6.0  DATA  REDUCTION 

4 


Several  different  data  reduction  methods  were  tried  in  searching  for  the 
best.  The  first  method  was  the  same  as  that  used  by  Cooke.  Cooke  recorded  data 
in  a  fashion  similar  to  the  procedure  mentioned  previously  and  used  Westphal's 
data  reduction  codes,  modified  to  account  for  the  pressure  gradient,  to  reduce  the 
data  he  obtained.  Although  the  error  criterion  determined  by  Westphal  produced 
results  within  ±5%  uncertainty  only  for  his  simpler  flow,  his  criteria  were  also 
applied  to  this  flow. 

The  second  method  used  the  derivative  of  the  oil  film  height  with  respect 
to  time  to  determine  the  skin  friction  since  the  time  rate  of  change  of  the  oil-film 
height  is  the  directly  measured  quantity.  Rearranging  equation  (1)  to  isolate  the 
height 


h 


r  t 


(31) 


and  taking 


4 


the  derivative  of  the  oil  film  height  with  respect  to  time  yields 
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31 

Since  is  inherently  negative  for  an  oil  film  thinning  under  the  action 
ot 

of  a  shear  stress,  multiplying  through  by  -  I  and  taking  the  logarithm  of  each  side 
yields  a  linear  equation 


log(  -  -f”  )  =  lQg(  -T-  )  “  2  io8(') 

y  —  L  +  mx 


(33) 

(34) 


r)h 

By  evaluating  -7-  as  a  backward  finite  difference  and  substituting 
ot 

h  -  hxN ,  discrete  values  can  be  determined  for  log^  -  J  and  log(f).  Thus  the 
equation  is  in  the  form  of  a  straight  line  with  slope  -2.  Using  a  linear  least 
squares  fit,  the  constant  term,  log^  ~~  )  =  b  ,  and  the  related  uncertainties  for 
each  beam  were  determined.  The  slopes  for  most  data  records  were  between  -2 
and  -2.2.  The  errors  on  the  constant  terms  were  about  3%  on  average.  Solving 
for  the  skin  friction  and  eliminating  x0  as  before,  the  final  equation  becomes 


[x2  —  10(fel-&2>x1] 
10M  -  I0bl 


(35) 


The  values  of  the  constants  b\  and  b2  were  fairly  close,  thus  th  ir  position  in  the 
denominator  makes  the  skin  friction  very  sensitive  to  these  values.  In  addition, 
by  taking  10  to  the  power  of  b\  and  bl  ,  any  errors  associated  with  these 
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constants  are  increased  accordingly  and  degrade  the  usefulness  of  the  calculated 
skin  friction.  This  process  results  in  an  unacceptably  high  uncertainty  of  as  much 
as  50%  on  the  final  result. 

The  third  method  was  actually  a  modification  of  the  first  method.  It  was 
realized  that  more  information  could  be  obtained  from  each  data  record.  Most 
notably,  the  negatively  sloped  zero  crossings  were  used.  These  points  were  not 
used  in  the  calculation  of  C(.x)  and  therefore  can  be  used  to  obtain  another  value 
of  C(jc)  and  shear  stress  using  data  from  the  same  test  run.  In  addition,  the  data 
was  shifted  up  by  .1  rms  of  each  signal  and  the  positive  and  negative  zero 
crossings  were  used  to  calculate  two  more  values  of  the  shear  stress.  The  data 
was  then  shifted  down  from  the  original  level  by  .1  rms  of  each  signal  and  the 
process  was  repeated.  Thus,  it  was  possible  to  obtain  six  values  of  the  skin 
friction  from  a  single  test  run,  which  translated  into  18  possible  independent 
values  of  skin  friction  at  every  location,  in  each  of  the  three  directions. 

Applying  the  third  method  produced  large  amounts  of  scatter  among  the 
skin  friction  values.  In  fact,  measurements  taken  at  the  same  location  and  in  the 
same  direction  varied  so  much  that  any  effort  to  combine  them  with 
measurements  in  other  directions  at  the  same  location  sometimes  resulted  in 
uncertainties  on  the  final  values  in  excess  of  100%. 

The  fourth  and  final  method  extended  the  third  method,  using  another 
piece  of  physical  information,  to  reduce  the  scatter  among  the  skin  friction  values 
calculated  at  a  given  point  in  a  given  direction.  In  reviewing  the  output  from  the 
reduction  program  and  seeking  a  method  of  reducing  the  uncertainty  of  the 
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results,  a  pattern  between  the  magnitude  of  the  skin  friction  and  the  time  origin 
error,  dt,  within  each  run  was  ascertained.  Many  plots  were  made  of  skin  friction 
magnitude  versus  dt  at  a  given  location  and  angle.  Figures  26,  27,  and  28  are 
typical  of  these  plots.  The  data  points  from  each  individual  run  are  connected 
with  straight  lines.  For  the  large  majority  of  test  runs,  the  magnitude  of  the  skin 
friction  component  is  almost  a  linear  function  of  dt  within  each  run.  Different 
test  runs  at  the  same  angle  do  not  Fit  onto  the  same  line,  or  even  necessarily  have 
slopes  of  the  same  sign.  However,  in  well  over  half  of  the  locations  viewed,  if  the 
lines  formed  by  connecting  the  values  from  each  test  run  were  extended  to 
dt  =  0,  the  values  converged. 

It  was  realized  that,  since  the  time  origin  is  a  property  of  the  particular 
flow,  and  therefore  independent  of  xf  or  x2,  the  time  origin  for  each  beam  should 
be  identical.  Westphal's  program  included  the  option  to  either  input  the  time 
origin  for  each  beam,  or  to  allow  the  program  to  calculate  the  time  origin  for  each 
beam.  By  varying  the  number  of  zeroes  used  from  each  beam,  it  was  normally 
possible  to  obtain  a  time  origin  error  within  the  stated  criteria  {dt  <  .02).  In  fact, 
in  most  cases,  it  was  possible  to  obtain  dt  values  less  than  .01.  Once  this  had 
been  accomplished  by  interactively  iterating  on  the  number  of  zeroes  for  each 
beam,  an  average  of  the  upstream  and  downstream  beam  time  origins  was  input 
as  the  time  origin  for  each  beam.  Thus,  the  time  origin  was  forced  to  be  the  same 
for  each  beam  and  the  time  origin  error  was  identically  zero.  This  process  usually 
required  between  3  and  5  iterations.  If  the  errors  on  the  constants  C(xt)  and 
C(x2)  were  still  below  the  acceptable  level  (e<.004),  then  the  skin  friction  values 
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computed  were  used  in  determining  the  magnitude  and  direction  of  the  skin 
friction  at  that  location.  A  typical  example  of  this  iteration  process  is  shown  in 
Table  1. 

The  first  column  of  Table  1  gives  the  number  of  zeroes  used  from  the 
upstream  beam  /  downstream  beam.  The  second  column  shows  the  skin-friction 
coefficients  obtained  in  the  first  reduction  iteration.  The  top  line,  denoted  "0  =", 
is  the  reduction  using  the  positive  sloped  zero  crossings.  The  second  line,  denoted 
"M  =  ",  is  the  value  obtained  by  subtracting  .1  of  the  signal  rms  from  the  average 
signal  intensity  and  using  the  positive  sloped  zero  crossings.  The  third  line, 
denoted  "P  = ",  is  the  value  obtained  by  adding  .1  of  the  signal  rms  value  to  the 
signal  average  and  using  the  positive  sloped  zero  crossings.  In  the  next  three 
lines,  the  first  letter  again  denotes  the  value  used  for  the  signal  mean  and 
calculated  the  skin  friction  using  the  negative  sloped  zero  crossings. 

In  Table  1,  part  (a),  the  data  records  (upstream  and  downstream  beams) 
were  reduced  using  20  zero  crossings  from  each  beam  and  the  time  origins  were 
calculated  by  the  computer  to  obtain  a  best  fit  to  equation  (14).  Note  that  the 
lowest  value  of  the  skin  friction  coefficient  was  .00175  and  the  highest  value  was 
.00187.  Also,  some  of  the  time  origin  errors  ( dt )  were  above  the  acceptable  limit 
of  .02.  Therefore,  the  data  were  reduced  again  using  10  zero  crossings  on  the 
upstream  beam.  The  output  is  shown  in  Table  1,  part  (b).  Again,  the  computer 
was  allowed  to  calculate  individual  time  origins  for  each  skin  friction  coefficient. 
Reducing  the  data  using  these  parameters  produced  time  origin  errors  within  the 
acceptable  limits.  Therefore,  10  zero  crossings  were  used  on  the  upstream  beam 
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and  20  zero  crossings  were  used  on  the  downstream  beam.  The  time  origin  was 
fixed  at  75.5  (the  average  value)  for  each  beam  for  the  final  run.  Using  these 
parameters  produced  the  output  in  part  (c).  Notice  that  the  maximum  difference 
in  the  skin  friction  coefficient  values  for  the  final  run  is  much  less  than  1%  of  the 
skin  friction  coefficient.  This  method  increases  the  precision  of  the  skin  friction 
coefficient,  but  not  necessarily  the  accuracy.  It  is  possible  that  another 
combination  of  the  number  of  zero  crossings  used  from  each  beam  could  also 
produce  an  acceptable,  but  different,  skin  friction  coefficient. 

Once  all  of  the  individual  data  records  had  been  reduced  with  dt  —  0  , 
factors  such  as  the  uncertainty  in  the  angle  between  the  oil-film  leading  edge  and 
the  measurement  direction,  and  the  effects  of  non-zero  shear  gradients  had  to  be 
accounted  for. 

The  analysis  at  the  end  of  Chapter  3,  in  which  equation  (29)  was 
developed  for  unknown  flow  angles,  was  based  on  the  perpendicularity  of  the 
initial  oil  line  and  the  measurement  direction.  Since  the  oil  line  was  applied  to  the 
surface  by  hand,  this  angle  could  not  be  precisely  ensured.  The  effects  of  errors 
in  angular  position  may  be  determined  through  an  analysis  of  Figure  29.  This 
figure  depicts  an  oil  flow  with  an  unknown  direction  y  with  respect  to  the 
measurement  direction  and  an  oil  line  which  deviates  from  being  perpendicular 
to  the  measurement  direction  by  an  angle  /J.  The  distance  travelled  by  an  oil 
particle  from  the  leading  edge  to  point  x,  is  —  5b,  and  the  distance  to  x2  is 
s2  —  Thus,  the  basic  shear  stress  equations  are 
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(36) 


,  _  Pifi  -Sqi)  ^  -  ^o)  cos  y  +  (x,  -  Jy0)  sin  y  tan(y  -  /?)] 

h[ti  hltl 

_  m(^2  -  *02)  _  v[(x2  -  *o) cos  y  +  (*2  -  *o)  s»n  V  tan(y  -  /?)] 
h2t2  h2t2 


Eliminating  xQ  as  before  and  reducing  yields 

r  i  i 

h2t2  -  h[tl  1  +  tan  y  tan  /J  cos  V 


(38) 


and 


t'x  =  r'  cos  y  =  -r* 


1 

1  +  tan  y  tan  /? 


(39) 


where  rx  is  the  shear  stress  component  in  the  measurement  direction  calculated 
from  equation  (17)  and  t'  is  the  total  shear  stress  (with  correction  for  angle  (l ). 

However,  /?  is  an  unknown  error  angle  so  the  total  shear  stress  cannot  be 
computed.  The  effect  of  this  error  factor  increases  as  the  angle  between  the 
measurement  direction  and  the  flow  angle,  y  increases.  When  the  measurement 
angle  and  flow  angle  coincide  (y  =  0)  the  error  factor  is  1.0,  and  therefore  the 
angle  P  has  no  effect.  Since  this  effect  was  not  considered  until  after  the  data  had 
been  obtained,  the  uncertainty  in  the  angle  P  was  not  minimized.  However, 
under  normal  conditions  (i.e.  sitting  at  a  desk),  the  error  in  this  angle  is  less  than 
±  1°  if  care  is  taken.  However,  for  purposes  of  data  reduction,  this  error  was 
assumed  to  be  less  than  ±  3°  to  account  for  the  limited  working  space  in  the  wind 
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tunnel  while  applying  the  oil  (see  Figure  14).  The  uncertainty  due  to  this  angle 
was  included  in  the  data  reduction  program  LS.FOR. 


The  geometry  for  the  shear  gradient  analysis  (see  Figure  7)  is  similar  to 
that  for  the  misalignment  angle,  p,  analysis.  In  all  previous  analyses,  the  shear 
stress  was  assumed  to  be  constant  and  was  denoted  simply  as  t.  Since  the  shear 
stress  is  now  assumed  to  be  variable,  t,  and  t2  are  used  to  represent  the  shear 
stresses  at  points  ($,,«,)  and  (52,  n2),  respectively. 

HO) 

"1*1 


T  •>  = 


(*[s2  ~  502) 


hit 


2*2 


(41) 


Furthermore,  expanding  the  shear  stress  in  a  Taylor  series  about  point  (s,n,) 
yields 


T’  =  T',  +  -jgr  (x'  -x',)  +  -~r  (z'  -  Z',)  + ... 


(42) 


Keeping  only  the  first  order  terms,  writing  the  above  equation  for  r{s^n 2)  =  r2  and 
substituting  for  ?!  and  t2  from  equations  (40)  and  (41)  yields 


~  50l)  ~~  5oi)  dr  (,  .  x  ,  dr  ,  ,  ,  , 

+  1F(  2“  l)  +  _a7(  2_ 


h-,t 


2*2 


hit 


1*1 


(43) 


Making  the  following  geometrical  substitutions 
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*1  -*0 

51  %  cos  y 


(44) 


*2  “  *0 

s2  ^02  COS  y 

(45) 

x>2  ~  x>  l  =  (*2-*l)C0S  CLd 

(46) 

z'2  -  z\  =  (x2  -  Xj)  sin  ad 

(47) 

Solving  for  x0,  substituting  into  equations  (40)  and  (41),  and  multiplying  each  by 
cos  y  to  obtain  the  shear  stress  component  in  the  direction  of  measurement  x 
yields 


~  *i) 

h2t2  -  Vi 


(48) 


r-y-,  = 


h2t2  -  hltl 


(49) 


However,  the  laser  beams  were  positioned  such  that  the  data  point  was  midway 
between  the  two  beams.  Thus,  since  only  the  linear  terms  were  retained  in  the 
Taylor  series  expansion,  the  value  of  the  shear  stress  at  the  data  point  is  taken 
as  the  average  of  r,  and  t2 


*  *  =  T* 


1  - 


h\tx+h2t2  cosy 


^cos^+-^sin“‘' 


(50) 
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Rearranging  equations  (15)  and  (16)  as  follows 


^1  -  hA^effl  (effl  -  \ — —  (5 1 ) 

hlh  =  ^A^eff2{efp.  =  ~~2  x - -  (52) 

shows  that  the  constants  /?,r,  and  h2t2  are  proportional  to  the  distance  from  each 
of  the  beams  to  the  leading  edge  of  the  oil  film.  Thus,  the  effects  of  non-zero 
shear  gradients  can  be  minimized  by  keeping  the  distance  from  the  upstream 
beam  to  the  leading  edge  (jc1  —  x0)  as  small  as  possible.  This  minimization, 
however,  is  limited  by  the  fact  that  x,  must  be  far  enough  away  from  x0  to  obtain 
sufficient  fringes  for  data  reduction.  Westphal  (1986)  recommends  using  at  least 
10  fringes  for  the  upstream  beam  calculations.  Decreasing  the  distance  between 
the  beams  ( x2  —  x,)  also  reduces  the  distance  between  the  downstream  beam  and 
the  oil  film  leading  edge,  which  reduces  the  effect  of  non-zero  shear  gradients. 
Therefore,  the  minimum  laser  beam  separation  required  to  maintain  spatial 
resolution  should  be  used. 

If  the  locations  at  which  data  were  taken  composed  an  orthogonal  grid,  the 

shear  gradients  could  be  obtained  through  an  iterative  procedure.  However,  since 

this  was  not  the  case,  a  method  was  needed  to  determine  the  shear  gradients  with 

the  given  data.  Noting  that  ( rr-r  cos  ad  +  -77  sin  ad)  from  equation  (50) 

ox  oz 

represents  the  shear  gradient  in  the  direction  of  measurement  (~^~)  »  an 
approximation  of  the  shear  gradient  can  be  obtained  from  the  data  records. 
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Referring  to  the  linear  oil-film  profile  in  Figure  4,  the  following  result  may  be 
obtained  via  simiiar  triangies  at  time  r, 


_  Mfi)  _  Mfi) 

x2  -xl  x2  ~  x0  jc,  -  x0 


(53) 


where  h,  and  h2  are  the  heights  of  the  oil  film  at  x,  and  *2  at  the  same  time.  The 
easiest  time  to  determine  the  height  of  the  oil  film  is  at  one  of  the  zero  crossings, 
but  the  upstream  and  downstream  beams  do  not  cross  the  zero  axis  at  the  same 
times.  Therefore,  the  height  of  the  upstream  beam  was  determined  at  one  of  its 
zero  crossings,  and  the  height  of  the  downstream  beam  was  found  by  linear 
interpolation  at  the  same  time.  Once  h,  and  h2  are  known,  the  distances  x,  -  x0 
and  x2  -  x0  can  be  determined  from  equation  (53).  Thus,  the  shear  at  the 
upstream  and  downstream  locations  can  be  determined  as  follows 


_  ~  *o) 

Vi 


f*(x2  ~  xo) 

T2  = -  - 

h2t2 


The  shear  gradient  was  then  taken  as 


dr  _  t2  Tl 
dx  ~  x2  ~  xl 


(54) 


(55) 


(56) 


Having  reduced  the  individual  data  records  and  accounted  for  uncertainty 
on  these  values,  the  skin  friction  values  that  met  the  error  criteria  were  combined 
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to  give  the  magnitude  and  direction  using  program  LS.FOR  (Appendix  C). 
Program  LS  comoines  the  acceptable  skin  friction  measurements  in  a  least  square 
error  fit  of  the  data  to  obtain  the  x'  and  z'  components  of  the  skin  friction 
coefficient.  Each  of  these  skin  friction  coefficients  should  satisfy  the  following 
equation 


~  Cx'  +  C2>  sin  adj 


(57) 


and  Q  were  subsequently  chosen  to  minimize  the  square  of  the  errors 


5  =  y  {Cmjl  -  Cx>  cos  ad4  -  C2>  sin  udj)2  =  minimun > 


4=1 


(58) 


To  obtain  the  values  of  Cy  and  Ce  wnich  yield  a  minimum  for  S,  the  partial 
derivatives  of  5  with  respect  to  Cy  and  Ct  were  set  equal  to  zero 


dS 

dCx. 


=  0  = 


9 

cos  ctdj  +  Cx>  cos2adt/  +  C2>  cos  <xdd  sin  adj) 

i=i 


(59) 


dS 

8Cy 


9 

=  0  =  ^(-cm,  sin  <xdd  +  Cx>  cos  add  sin  a dd  +  C2>  sin2a^) 

4=1 


(60) 
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With  the  following  definitions 


9 

-Z 


A  =  ZjC0S  “ dJ 

i=l 


(61) 


y 

B  =  £cos  udd  sin  adJ 


1=1 


(62) 


9 

-z 

1=1 


sin  ctdJ 


(63) 


T 

£  =  cos 


1=1 


(64) 


F=  ^Cmj  Sin 


i=i 


(65) 


these  equations  may  be  written  in  matrix  notation 


‘ A  B 1 
.B  C 


c  1 


'£l 


(66) 


The  values  of  Q  and  C,  which  best  fit  the  data  Q,  are  then  given  by 
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DE-BF 


(67) 


Q  = 


AD-  B' 


<=Z'  = 


AF- BE 
AD  -  B2 


(68) 


These  components  were  then  combined  to  give  the  magnitude  and  direction  of  the 
skin  friction  coefficient 


Cf-yfe  +  C} 


(69) 


a  =  tan 


l  Cx> 


(70) 


6.1  Uncertainty  Analysis 


Program  LS  averaged  all  six  values  from  a  single  run  together  and  computed  95% 
uncertainty  bounds  on  the  averaged  value  from  each  run.  Thus,  for  the  general 
case  in  which  there  were  three  test  runs  taken  in  each  of  the  three  measurement 
directions,  LS  first  computed  nine  skin  friction  coefficients  along  with  their 
uncertainties.  Then  the  uncertainty  in  each  measurement  due  to  the  leading  edge 
perturbation  angle  /J  was  added  to  the  original  uncertainty  by  taking  the  square 
root  of  the  sum  of  the  squares.  This  produced  nine  skin  friction  coefficients  and 
their  related  uncertainties  --  Cra  and  <5 Cm,  i  =  1,2,. ..,9.  The  least  squares  fit  was 
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then  performed  as  described  above  to  obtain  and  Q.  These  values  were 
considered  the  "baseline"  values. 

To  determine  the  uncertainty  in  C*  and  Q  due  to  each  of  the  individual 
measurements,  the  first  value  was  perturbed  by  an  amount  equal  to  its 
uncertainty 


C'm.i=CmA+dCmA  (71) 

and  the  least  square  error  fit  was  performed  again.  The  differences  between  the 
Ct  and  Ct  obtained  from  the  perturbed  calculation  and  the  baseline  case  were 
taken  as  the  uncertainties  in  and  Q  due  to  the  uncertainty  in  the  first 
measured  skin  friction  coefficient  CmA.  Each  of  the  remaining  eight  values  were 
similarly  perturbed  to  obtain  the  uncertainties  in  C,  and  C,  due  to  the 
uncertainty  in  each  of  the  original  values  Cm>,  These  nine  uncertainties  were  then 
combined  with  the  least  square  error  S,  by  taking  the  square  root  of  the  sum  of 
the  squares  to  produce  the  final  uncertainties  on  Cy  and  Q  of  <5Q  and  <5Q, 
respectively.  For  comparisons  in  Table  2,  the  magnitude  and  angle  were 
necessary,  and  thus  the  uncertainties  needed  to  be  expressed  in  terms  of  the 
magnitude  and  angle.  Therefore,  bCt  and  <5Q  were  combined  as  follows  to 
obtain  bCf  and  b a 

^  <72) 


where 
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(73) 


from  equation  (69),  and 


from  equation  (69) 


where 


from  equation  (70),  and 


cl  +  cl 


(77) 


from  equation  (70). 
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6.2  Curve  fitting  method 


An  alternative  method  for  reducing  the  data  was  developed  based  on  curve 
fitting  the  entire  data  record.  This  method  uses  every  element  in  the  data  record, 
instead  of  only  using  the  zero-crossings  or  signal  peaks,  to  obtain  the  skin  friction. 
Knowing  that  the  intensity  of  the  returning  laser  beams  varies  sinusoidally  as  a 
function  of  the  oil  film  height  and  that  a  change  in  height  of  hk  corresponds  to 
one  period,  the  intensity  may  be  written  as 


A  sin 


(78) 


where  A  is  the  amplitude  of  the  signal.  As  r -*  oo,  h-*0,  and  1.  Thus, 

/ 1 

<f)  =  -j-  and  the  expression  for  the  intensity  becomes 


A 


(79) 


Solving  equation  (4)  for  h  and  substituting  yields 

where 


k  = 


n(x  -  -x0) 

T  t 


(80) 


(81) 
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The  constant  k  cannot  be  found  by  simply  solving  equation  (80)  for  k  and 
evaluating  the  expression  for  all  -j-  since  cos-1^-^-^  oscillates  between  -1  and  1 
while  t  increases  continually.  Multiple  values  of  k  yield  the  same  value  of  in 

r I 

equation  (80)  since  cosine  is  a  multiple-valued  inverse  function.  Thus,  it  is 
necessary  to  add  2n  to  the  value  of  k  for  each  successive  period.  In  addition, 
unless  the  oil  is  allowed  to  thin  to  less  than  one  fringe  in  height,  the  multiple  of 
2 n  to  be  added  is  unknown.  After  careful  inspection  of  this  phenomenon,  the 
following  equation  was  derived  which  can  be  iteratively  solved  for  k 


k  =  ~2xS  ‘os-'(AW, 


where 


+  lfor  m  <  <  m  + 

—  1  for  m  +  <  -y-  <  m  +  1 


P  —  integer  portion  of 


The  value  of  S  may  be  represented  analytically  by 
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5  = 


(85) 


The  iteration  on  equation  (82)  is  a  two  variable  iteration  (on  k  and  t0  since  the 
time  t  in  equation  (82)  is  referenced  to  the  time  origin  /0 

However,  the  current  method  uses  only  the  data  elements  at  the  zero 
crossings  to  determine  the  skin  friction  component  from  each  data  record.  This 
results  in  using  only  20  or  30  data  elements  of  the  4000  recorded.  Since  the  curve 
fitting  procedure  uses  every  data  element  recorded,  it  is  likely  to  reduce  the  total 
number  of  data  elements  necessary  and  thereby  reduce  the  amount  of  time 
necessary  to  record  data. 
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7.0  DISCUSSION 


The  locations  at  which  data  were  obtained  fall  into  six  basic  groups  plus 
a  few  miscellaneous  points.  Four  of  the  groups,  consisting  of  points  10-17,  20-27, 
31-38,  and  41-49,  were  arranged  such  that  all  the  points  within  a  group  are 
located  on  a  line  extending  away  from  the  wing-body.  The  lowest  numbered  data 
point  in  each  of  th-0e  groups  was  furthest  away  from  the  wing  body,  and  the  rest 
count  inward  toward  the  wing- body.  Points  1-9  and  61-66  were  located  in  the 
nose  region  and  are  all  located  outside  the  line  of  low  shear.  Points  51  and  52 
were  located  far  from  the  wing-body  along  the  side.  Most  of  the  remaining  data 
points  were  located  in  the  plane  of  symmetry  of  the  wing-body. 

Skin  friction  coefficient  vectors  measured  at  the  various  positions  in  the 
wing-body  junction  flow  are  plotted  in  Figures  30  and  31,  and  are  presented  in 
tabular  form  in  Table  2.  The  first  three  columns  of  Table  2  contain  data  location 
information  --  identification  number,  and  the  x ',  z'  coordinates 
non-dimensionali2ed  by  the  maximum  wing  thickness  T.  The  next  two  columns, 
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labeled  Cf  and  a,  contain  the  magnitude  and  direction  of  the  skin-friction 
coefficients  obtained  using  the  oil  film  interferometry  method  discussed  in  the 
previous  chapter.  The  initial  height  of  the  oil  films  (h  =  h^N0i2)  used  in  this  study 
were  between  10 nm  and  20 fim. 

The  column  labeled  aMJ  contains  flow  directions  which  were  measured 
directly  from  77C>2-K.erosene  surface  oil  flow  visualizations  similar  to  that 
represented  in  Figure  2.  To  obtain  these  angles,  a  tangent  was  drawn  to  the  oil 
streaks  at  the  point  of  interest,  and  the  angle  of  the  tangent  with  respect  to  the 
tunnel  centerline  was  measured.  These  angles  could  only  be  determined  to  within 
±  3°  in  most  locations  with  significantly  higher  uncertainties  (±  10°)  in  the  nose 
region  (points  13-17).  The  angles  measured  from  the  surface  oil  flow  visualization 
agree  to  within  the  experimental  error  of  the  angles  measured  using  the  Oil- Film 
Laser  Interferometry  technique  in  other  areas  away  from  the  line  of  low  shear. 
Near  the  line  of  low  shear,  the  angles  show  the  greatest  disagreement.  The  heights 
of  the  oil  films  used  in  these  flows  were  less  than  20  nm.  However,  the  oil  used 
in  the  surface  oil  flow  visualization  was  several  times  thicker,  and  varied  in 
thickness  inversely  with  the  shear.  Thus,  these  angles  are  more  likely  to  be 
affected  by  pressure  gradients. 

The  column  labeled  C/IO|  contains  skin  friction  values  at  most  of  the 
locations  examined  in  this  thesis.  These  skin  friction  values  were  obtained  by 
fitting  the  semi-logarithmic  region  of  mean  velocity  profiles  (measured  using  a 
laser  anemometer)  to  the  log  law  for  equilibrium  boundary  layers  (Devenport  and 
Simpson,  1989)  according  to  the  equation 
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.  yu,  ITT  +  w2 

where  y  +  =  — — ,  wt  =  yj  —  ,  and  u  +  was  taken  to  be - - - .  However, 

since  the  log  law  is  a  two-dimensional  formulation,  this  method  was  not  expected 
to  yield  reliable  values  of  Cf  under  the  vortex  legs,  but  provide  a  qualitative  basis 
for  comparison. 

The  skin  friction  coefficients  obtained  by  fitting  the  semi-log  region  and 
those  obtained  by  oil  film  interferometry  are  plotted  as  a  function  of  ID  number 
for  comparison  in  Figure  32.  in  the  region  in  front  of  the  nose,  both  sets  of  skin 
friction  coefficients  exhibited  the  same  trend,  although  the  C/log  values  exhibit 
much  greater  fluctuations.  Near  the  maximum  wing  thickness,  the  trends  are  also 
alike,  C,  increases  as  the  wing  is  approached,  with  a  drop  in  the  magnitude  near 
the  line  of  low  shear  (located  in  the  proximity  of  points  23  and  24).  In  this  area, 
however,  the  values  of  C/log  are  much  higher  than  those  measured  using  oil  film 
interferometry.  In  the  regions  along  the  side  and  rear  of  the  wing,  the  trends  are 
again  similar,  and  the  discrepancy  between  the  magnitudes  decreases  with 
distance  downstream.  Notice  that  these  discrepancies  between  Cf  and  C/log 
increase  rapidly  inside  the  line  of  separation.  This  trend  is  expected  since  the  log 
law  is  not  expected  to  hold  under  the  vortex  leg.  Furthermore,  it  the  differences 
between  the  two  methods  is  smaller  downstream,  where  the  vortex  strength  has 
decreased,  than  upstream. 
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The  close  agreement  away  from  the  vortex,  adds  credence  to  both  methods  for 
such  a  flow.  However,  since  the  C/log  values  can  only  be  considered  as  a 
qualitative  guide  under  the  vortex,  the  similar  trends  exhibited  by  the  Cf  and  the 
C/log  methods  is  encouraging.  Sources  of  error  due  to  three-dimensional  effects 
will  be  discussed  later  in  this  chapter. 

The  column  labeled  aw/u  contains  flow  angles  obtained  by  combining  x'  and 
z'  velocity  components  which  were  measured  by  hot-wire  anemometry  (points 
1-9)  and  laser  Doppler  velocimetry  (points  10-27  and  61-66).  The  angles  for 
points  10-27  were  computed  from  the  following  equation 

aw«  =  tan-'(-£-)  (87) 

The  angles  for  points  61-66  were  measured  by  Oicmen  (1990)  in  a  similar 
manner.  Velocity  component  data  were  not  available  at  other  locations  in  the 
flow.  The  height  at  which  the  velocity  components  were  measured  at  each 
location  are  presented  in  the  final  column,  labeled  hwlll.  The  angles  given  at 
locations  1  and  2  are  somewhere  in  the  peak  between  the  two  legs  of  a  Johnston 
(1960)  hodograph  plot.  Therefore,  the  actual  angles  at  these  locations  are  at  least 
as  large  as  those  given  in  Table  2,  and  probably  slightly  larger.  The  rest  of  the 
awlu  data  were  located  on  the  near-wall  leg  of  the  Hodograph  plot  and  can 
therefore  be  considered  accurate  as  near  the  wall  as  is  known.  Uncertainties  were 
available  only  for  points  61-66,  and  tend  to  be  very  large  near  the  wall 
(characteristic  of  LDV  angle  measurements  near  the  wall).  At  points  61-66,  the 
LDV  angles  and  the  oil-film  interferometry  angles  agreed  to  within  the 
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experimental  uncertainties  given.  The  velocity-component  angles  and  the  laser 
interferometry  angles  follow  the  same  general  Umd  -  increasing  as  the  flow 
approaches  the  wing,  attaining  a  maximum  near  points  10-17,  and  decreasing 
around  the  side  of  the  wing,  as  expected.  However,  the  oil  film  interferometry 
angles  exhibited  much  larger  variations  than  the  velocity  component  angles. 
However,  as  noted  previously,  as  the  wall  is  approached,  the  uncertainties  in  the 
LDV  angles  becomes  very  large  —  even  at  15 0  fxm  .  No  angles  were  available  for 
comparison  with  the  oil-film  interferometry  angles  which  were  measured  at  the 
same  height. 

At  several  points  in  the  flow  (marked  with  a  D  in  Table  1),  data  could  not 
be  obtained.  A  separation  point  exists  along  the  centerline  at  x'/T  =  —  .47. 
Points  29  and  30  are  between  this  point  and  the  model  nose.  In  this  region,  the 
direction  of  flow  is  from  the  nose  towards  the  separation  point,  and  the  flow 
direction  changes  rapidly  away  from  the  centerline  (Figure  33).  As  a  result,  the 
oil  does  not  flow  over  both  of  the  laser  beams.  The  situation  could  be  remedied 
by  significantly  decreasing  the  beam  separation.  Points  39  and  40  are  located 
behind  the  tail  of  the  model.  In  this  region,  the  oil  was  observed  to  diverge  from 
the  tunnel  centerline  (Figure  34).  This  divergence  likewise  prevented  data  from 
being  obtained  in  this  region  since  the  oil  did  not  flow  over  both  of  the  laser 
beams. 

At  points  marked  with  an  A,  B,  or  C  in  Table  1,  data  could  not  be 
obtained  in  one  of  the  three  directions  of  measurement.  This  occurred  when  the 
angle  of  measurement  was  nearly  normal  to  the  measurement  direction,  as 
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depicted  in  Figure  8.  This  data  loss  was  not  critical,  but  the  uncertainties  quoted 
are  questionable  since  the  largest  source  of  uncertainty  was  normally  the  error 
of  resulting  from  the  least  squares  Fit. 

Three  dimensional  effects  are  illustrated  in  Figure  30.  Note  that 
downstream  of  the  maximum  wing  thickness  the  flow  speed  and  direction  vary 
gradually.  The  exception  to  this  is  the  data  point  furthest  downstream  and 
closest  to  the  tunnel  centerline  (point  33).  At  this  location,  the  flow  diverges 
sharply  from  the  centerline  creating  large  shear  gradients.  At  points  upstream 
of  and  near  the  maximum  wing  thickness  three-dimensional  effects  are  larger  due 
to  the  flow  acceleration  around  the  wing  and  are  very  significant  close  to  the  wing 
(points  13-17).  effects.  The  flow  angles  change  significantly  over  the  distance 
between  the  laser  beams,  and  the  streamline  radii  of  curvature  are  on  the  order 
of  the  distance  between  the  laser  beams.  Since  the  laser  beams  were  positioned 
such  that  the  measurement  location  was  located  midway  between  the  laser  beams, 
the  laser  beams  impinged  on  the  test  surface  at  different  locations  when  aligned 
in  each  of  the  three  measurement  directions.  The  flow  varied  significantly  over 
the  distance  between  the  laser  beams.  Therefore,  skin  friction  measurements  in 
the  different  directions  (0°,  -45°,  and  -60°  )  were  measuring  skin  friction  in 
significantly  different  flows.  Attempts  to  combine  these  data  produced  absurd 
results.  The  angles  predicted  from  the  oil  interferometer  reduction  method  for 
points  13-17  were  between  0°  and  25°  although  the  actual  oil  pattern  was 
observed  after  the  test  runs  to  be  approximately  50°  to  60°  . 
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Another  possible  source  of  error  is  the  variation  in  thickness  of  the  initial 
oil  line.  In  all  of  the  previous  analyses,  it  was  assumed  that  the  initial  oil  line  was 
of  uniform  thickness,  and  that  as  the  oil  thinned,  lines  parallel  to  the  oil  film 
leading  edge  were  also  of  uniform  thickness.  If  one  end  of  the  initial  oil  line  is 
thicker  than  the  other  end  (as  shown  in  Figure  35)  and  the  angle  y  between  the 
measurement  direction  and  the  flow  angle  is  significant,  the  initial  conditions 
along  the  oil  film  leading  edge  are  not  identical.  Therefore,  jt0  cannot  be  assumed 
identical  for  each  of  the  laser  beams  and  cannot  be  eliminated  from  the  oil 
lubrication  equations.  However,  during  and  after  test  runs  in  which  the  initial 
oil  line  was  observed  to  have  uneven  thickness,  fringes  of  constant  thickness, 
visible  in  white  light,  were  qualitatively  observed  to  be  very  close  to  parallel  to  the 
oil  film  leading  edge  by  the  time  the  oil  thinned  enough  to  begin  obtaining  data. 
Thus,  this  non-uniform  oil  thickness  was  observed  to  be  a  transient  effect  which 
subsides  quickly. 

The  pressure  gradients  encountered  in  this  flow  were  found  to  have 
minimal  effects  on  the  oil-film  interferometry  method  used.  The  largest  pressure 
gradient  effects  occurred  in  the  region  near  points  13-17.  When  the  pressure 
gradient  correction  term  was  included  for  these  points,  the  skin  friction  value 
produced  from  a  given  test  run  differed  only  in  the  fifth  significant  digit. 
Therefore,  the  pressure  gradient  effects  were  not  included  in  the  reduction  of  the 
remaining  points. 

The  error  induced  by  the  unknown  leading  edge  angle,  /?  was  on  the  order 
of  3%  to  5%  on  the  individual  test  runs.  The  least  square  error  fitting  procedure 
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tended  to  minimize  this  effect  so  that  the  uncertainty  in  P  of  ±  3°  had  little  effect 
on  the  final  results.  The  largest  source  of  error  was  normally  the  disagreement 
between  measurements  in  the  three  different  directions  when  combined  in  the 
least  square  error  fit. 
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8.0  CONCLUSIONS 


In  this  study,  improvements  were  made  to  the  basic  dual  beam  oil  film 
laser  interferometer  system  which  increased  the  ease  of  applicability  and 
decreased  the  uncertainty  of  the  final  results  for  a  three-dimensional  flow.  This 
improved  system  was  used  to  produce  skin  friction  coefficients  around  a 
wing-body  (hull-appendage)  junction.  In  the  course  of  this  investigation, 
limitations  of  this  system  were  encountered  and  methods  for  minimizing  these 
limitations  were  discussed. 

The  laser  assembly  was  secured  in  one  of  three  pre-set  directions  atop  a 
movable  platform  which  in  turn,  was  mounted  on  a  perpendicular  set  of  rails. 
This  mounting  decreased  the  time  required  to  re-locate  the  laser  assembly  at  new 
data  point  from  45  minutes  to  less  than  5  minutes.  In  addition,  an  iterative 
procedure  for  accurately  determining  the  time  origin  for  each  oil  flow  was 
developed. 


CONCLUSIONS 


57 


The  current  interferometer  system  was  not  able  to  obtain  data  in  the 
centerline  of  the  wake  or  in  front  of  the  wing's  leading  edge  due  to  the  flow 
patterns.  In  addition,  this  method  produced  very  poor  results  (angles  off  by  more 
than  40°  at  some  locations)  in  the  presence  of  large  three-dimensional  effects. 

It  was  not  possible  to  make  direct  comparisons  to  vouch  for  the  validity 
of  these  results.  The  intent  of  this  study  was  not  to  confirm  other  results,  but  to 
produce  needed  skin  friction  values  in  a  complex  flow.  It  was  possible  to  obtain 
flow  angles  from  surface  oil  flow  visualizations,  hot  wire  data,  and  LDV  data  and 
compare  these  angles  to  the  angles  computed  for  the  skin  friction  coefficients. 
The  angles  computed  in  this  study  agreed  well  with  angles  from  the  other 
methods  in  areas  in  which  the  flow  was  nearly  two-dimensional.  In  areas  where 
large  three-dimensional  effects  were  present,  all  of  the  angles  generally  disagreed 
with  each  other  since  all  the  angles  were  measured  at  different  heights  in  the  flow 
and  the  flow  angle  varies  rapidly  near  the  wall  in  these  regions. 

Skin  friction  data  were  available  at  some  locations  in  the  flow,  was  not 
quantitatively  correct  under  the  vortex,  and  the  data  inside  the  line  of  separation 
could  only  be  considered  qualitatively.  Away  from  the  line  of  separation,  the 
values  produced  by  the  two  methods  agreed  to  within  experimental  uncertainty. 

Significant  improvements  in  the  current  system  could  be  obtained  by 
reducing  the  beam  spacing  ( x2  —  *0  and  the  distance  to  the  leading  edge  of  the 
oil  film  (x]  -  jc0).  These  distances,  to  a  large  extent,  determine  the  applicability 
of  the  method  to  the  flow.  If  the  flow  changes  significantly  over  these  distances, 
then  the  method  will  not  work.  Thus,  reducing  these  distances  will  increase  the 
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complexity  of  the  flow  that  can  be  studied  using  oil-film  laser  interferometry. 
Furthermore,  since  this  system  measures  the  average  shear  stress  over  the  area 
between  the  beams,  reducing  (x2  —  x{)  will  make  the  Final  result  approach  the 
actual  point  value  more  closely. 

The  author  has  several  general  suggestions  for  increasing  the  performance 
of  the  present  system  for  application  to  very  complex  flows,  such  as  the  one  in  the 
present  study.  Some  method  or  mechanical  device  to  increase  the  accuracy  of  the 
angle  of  the  initial  oil  line  will  reduce  the  error  due  to  this  effect.  Although  this 
effect  was  generally  very  small  compared  to  the  other  errors  involved,  in  order  to 
obtain  very  low  uncertainty  data,  this  effect  can  be  greatly  reduced  if  considered 
beforehand.  Also,  since  knowledge  of  the  shear  gradient  is  necessary,  taking 
measurements  on  an  orthogonal  grid  may  prove  useful.  Furthermore,  a  method 
to  make  fine  adjustments  of  the  movable  platform  would  be  helpful,  and  a  more 
solid  support  structure  than  the  wooden  frame  used  in  this  study  should  be 
constructed. 

In  further  studies,  closer  investigation  of  the  curve  fitting  method  of  data 
reduction  should  have  high  priority.  The  reduction  in  time  necessary  to  obtain 
enough  data  for  data  reduction  could  potentially  be  as  significant  as  the  reduction 
in  time  which  occurred  in  moving  from  a  single  beam  method  to  a  dual  beam 
method.  Also,  especially  important  if  the  curve  fitting  method  is  used,  the 
photodetectors  should  be  operated  in  the  photoconductive  mode  instead  of  the 
photovoltaic  mode  since  the  photovoltaic  mode  may  distort  the  shape  of  the 
interference  signal. 
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Table  1.  Sample  data  reduction  output 


Part  (a):  First  reduction  iteration 
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Part  (c):  Final  reduction  iteration  (dt  = 
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Tabic  2.  Skin  friction  coefficients  and  comparisons 
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-  No  data  in  45  °  direction 

-  No  data  available  (See  Discussion  Section) 


Tabic  2.  Skin  friction  coefficients  and  comparisons 


E 

a. 


o  o  o  o  o  o 

i/i  1^1 i/i  i/i  in 


„  0  o  d  d 
“  *+5+5+5+) 

'  "  ^  o  o  o  o 

°  o\  < NfSrlN 
Tt  ,  C4  ro  m  <N 


5r*  O  oo 

r  O  <N 

VJ  ©  ©  © 
O  O  O 


io*  on  — 

m  i  • 


k/*>  r<*) 

o  o 
©  o 


r-  rr  oc  o  *n  r** 

©  o  c  o  o  © 
o  ©  o  ©  o  © 


— rovcoo'O 


o  ©  ©  oo  —  O 


0  00  oo 

•  o 

*  •  4  °  *  .  •  . 

0  0  u-i  00  —  \C 

on  ; 

00  vo 

^o;--Nid 

— • 

d  re 

rtci",“<1((jN^ 

(N  — 

+1+5 

+>0  o 

-H+) 

+<*+5+5+1^*+, 

,,  ,,+5+5+5+5 

+,+lO  0  0  4 

°  O'  Cf 

•  0 

*  •  °  ’  °  .  «  . 

0  8  o»-m 

°*  —  vri 

®  Os 

^  d  ob  +  + 

—  —  f+i 
*  • 

—  ci 

O  -  Tf  c- ^  d  ^ 

Tf  OO  —  ^  (N 

N®  V® 

oN  oN  ox 
—  **)  O 
vS  rS 

C'-H+f+f 


IT 

d  O 
-H-H 


o'  o'  o 
ao  oo  voi  re 

-H+f-H-H 


v«v?v»v? 
o'  o'  o'  o' 
d  vO  oo  d 
d  —  —  re 

+5 +5+5 +5 


5?  S?  5?  5?  S?  S? 
C  vO  Cl  Cl  lO  f'. 
—  vri  re  c-i  re  ^ 

-H-H-H-H-H+l 


*r 

O' 

<N 

O’  VO 

ON 

•— * 

— . 

c*> 

<N 

©  o 

fO 

CM 

NO 

<N 

n 

© 

Tf 

IT  — 

^r 

<N 

o 

v-» 

00 

«B9 

oo  m 

ON 

O' 

r* 

fN 

V4"i 

KT\ 

m 

r-n 

O 

d  d 

<N 

m 

rn 

r-i 

rO 

"S’ 

rs  rn 

m 

m 

(N 

(N 

§ 

© 

© 

o  o 

s 

© 

© 

© 

© 

C 

©  © 

© 

© 

© 

© 

© 

© 

© 

© 

o  o 

o 

© 

© 

© 

O 

©  © 

© 

© 

© 

© 

© 

© 

o 

© 

© 

d  d 

© 

© 

© 

© 

© 

© 

©  © 

© 

© 

© 

© 

© 

© 

,  oooooooonnvi'nininoo 
t-ooooo  —  —  —  —  —  —  —  —  —  re«e 
‘■—c  OOOOOOOOOOOOOQd 
(HjOCOOOOOOOOOOOOCOOOOOOOOlOOO 
l  'Tcin-o,r<NO?'»i"C'n<Tinn 

ooooo  —  ~~dooodd  —  ~ 


o  o  o  o  o  o 

mcKNClO'C 
ah-CCO 
-T  vO  —  d  d 
(N-CffiCO 

———dod 


oooooooooooooooo 
Hooooooocooooooo-id 
*»«  —  —  —  —  —  dddddddddCT'ir 
X'C'O'O'C'Ccoaooooooooooooooo-oo 
-r-rrer-TTr  —  —  ~  —  —  —  —  —  —  dd 


©  ©  ©  ©  ®  ® 
o  -S  cr  Cl  3  rs 
d  d  —  d  ©  © 
CdffiOr'- 
r-  tt  o  —  cr  oo 


'rwwciciciociecirtcin-  oooooo 


„  a  q  c  «  <  es  SQQcee  cs 

^vOdoooo  —  dre^vnvodoooo  — 

gclCicld'TW'r^^,T,T'Tvr>0<n 


u 

-dci^m 
'C  <0  nC  ^  'C 


3 


e  e 

o  .2 


«  u 

.g.s 

13  13 


.S.S 

<9  <« 
•;  — 
<9  rs 

•o  *o 
c  O 
ZZ 


c 

o 

I 

e 

,o 

*55 

<n 

3 

a 
cS 
•2  „ 
u  u 

1  w 

13 

•  J 3 
V>:§ 

*»  a 

.S  « 

2  ~ 

-3-3 

0  o 
ZZ 


<BUQ 


64 


Illustrations 


Illustrations 


65 


mam 


/ 


).003l 

/ 


.0038 


0.2147 


Figure  4.  Two-dimensional  oil  flow  nomenclature 
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Figure  5.  Geometry  for  optical  pathlength  difference 
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Figure  9.  VPI&SU  Low  Speed  Boundary  Layer  Wind  Tunnel 
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Figure  11.  Wooden  frame  surrounding  the  wind  tunnel 
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Figure  16.  Laser  assembly  viewed  from  optical  rail  axis 
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Figure  17.  Laser  assembly  geometry  diagram 


Figure  18.  Box  diagram  of  hardware  components 
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Sample  photodetector  output:  Point  05,  Angle  60,  Take 
downstream  beam 
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•  imc  (minuter.) 

Figure  22.  Sample  photodetector  output:  Point  II  AhpIp  ds 

upstream  beam  '  An8lc  45>  Ta^e 


mple  "bad"  photodetector  output:  Point  28,  Angle  45,  Take 
upstream  beam 
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Figure  25.  Sample  'bad'  photodetector  output:  Point  28,  Angle  45,  Take 
1 ,  downstream  beam 
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igure  26.  Skin  friction  coefficient  magnitude  vs.  time  origin  error:  Point 
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Figure  27.  Skin  friction  coefficient  magnitude  vs.  time  origin  error:  Point 
04,  angle  60 
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Figure  28.  Skin  friction  coefficient  magnitude  vs.  time  origin  error:  Point 


Figure  29.  Non-perpendicularity  Geometry  and  Notation:  Error 
angular  position  of  oil  film  leading  edge 
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Skin  friction  coefficient  vectors  plotted:  entire  flowfield 
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Appendix  A.  PROGRAM  OILAD32.BAS 


1000  ' 

1005  DEFINT  A-Z 

1010  OIL-FILM  INTERFEROMETER  A/D  PROGRAM.  VERSION  3.2  082587 
1020  ' 

1030  A/D  PROGRAM  FOR  TWO  BEAMS  &  INPUT  CHN’LS,  W,0  SCANNER  &  EXT  TRIG 
1040  GOSUB  20000 

1050  INPUT  PARAMETERS  FOR  DATA  FILE  AND  A/D  CONVERSION 
1080  CLS:PRINT:PR1  NT"  OIL  A. D  PROGRAM  &  DATA  FILE  PA RAMTERS*: PRINT 
1120  PRINT’  A.D  PROGRAM  FOR  DUAL  BEAM  LASER  INTERFEROMETER  '  PRINT 
1130  PRINT'  FOR  USE  W/O  SCANNER  OR  EXT  TRIGGER  (VER.  3.1)  -.PRINT 
1140  PRINT:INPUT-  X-LOCATION  OF  MEASUREMENTS  -;X! 

1145  PRINT  INPUT'  Z-LOCATION  OF  MEASUREMENTS  '.ZI 
1150  PRINTINPUT'  REF  PRESSURE-Q  (IN.  WATER)  ";UREF! 

1155  PRINTINPUT'  NUMBER  OF  DATA  SAMPLES  '-.NUM.CONV 

1156  PRINTINPUT'  DATA  SAMPLING  RATE  IN  HERTZ  *;HERTZa 

1160  PRINT  INPUT'  DISC  FILE  FOR  RESULTS  -;DS 

1 180  PRINTPRINT'  Ensure  <  CAPS  LOCK  >  is  on  before  continuing" 

1 190  PR1NT:PRINT:PRINT:PRI NT" Press  any  key  to  continue' 

1200  AS  =  INKEYS  AS  =  "THEN  1200 
1300  CLOSE  LOPEN  DS  FOR  OUTPUT  AS  ai 
1310  PRINTal, 'X  =  -;X! 

1320  PRINTal  ,"Z  =  ";Z! 

1330  PRINTal, -Q  =  '.UREFI 

1400  PERIODa  =  (ia/(2o*HERTZa)),1000000a 

1450  TICKSs  =CINT(PERIODa/1.25«-32768a)-^  32768a 

1460  PERIODa  =TICKSa*  1.25a 

1500  START.CHANNEL=  1:END.CHANNEL  =  2 

1600  NUM.CONV  =  ((END. CHANNEL-START. CHANNEL  +  l)*NUM.CONV)»2 

2000  GOSUB  21000 

3060  PARM  =  (NUM.CONV-4) 

4000  FOR  ADDRESS  =  0  TO  PARM  STEP  4.TPRINT  =  PEEK(ADDRESS)+  PEEK(ADDRESS -*■  l)*256! 
4005  TPRINTB  =  PEEK(ADDRESS  +  2)  +  PEEK(ADDRESS  +  3)*256! 

4010  PRINTa I. TPRINT, TPRINTB 
4020  NEXT 
4030  END 
20000  ’ 

20010  '  Initialise  A/D  converter 
20020  • 

20030  BASE.ADDRESS  =  &H2EC 

20040  COMMAND. REGISTER  =  BASE.ADDRESS  +  1 
20050  STATUS. REGISTER  =  BASE.ADDRESS  +  l 
20060  DATA.REGISTER  =  BASE.ADDRESS 
20070  COM  MAN  D.  WAIT  =  &H4 
20080  WR1TEWAIT  =  &H2 
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20090  READ. WAIT  =  &H5 

20100  CSTOP  =  &HF 

20110  CCLEAR  =  &H  1 

20120  CERROR  =  &H 2 

20130  CCLOCK  =  &H3 

20140  CSAD  =  &HD 

20150  CRAD  =  &HE 

20160  EXT. CLOCK  =  &H40 

20170  EXT  .TRIGGER  =  &H80 

20180  CDMA  =  &H10 

20190  DUMMY  =  5 

20200  DMACHANNEL  =  1 

20210  DMAMODE  =  &H45 

20220  BASER  EG  =  2 

20230  COUNTREG  =  3 

20240  PAGEREG  =  &H83 

20250  STARTPAGE  =  &H3000 

20260  DEF  SEG  =  STARTPAGE 

20270  DMABASEL  =  &H0 

20280  DMABASEH  =  &H0 

20290  DMAPAGE  =  STARTPAGE/&H  1000 

20300  ' 

20310  Stop  and  clear  the  DT2801. 

20320  OUT  COMMAND.REGISTER,  CSTOP 
20330  TEMP  =  INP(DATA. REGISTER) 

20340  WAIT  STATUS. REGISTER.  COMMAND. WAIT 
20350  OUT  COMMAND.REGISTER,  CCLEAR 
20360  ' 

20370  CLS:PRI.NT:PRINT*  SETTING  UP  INTERFACE* 

20410  ' 

20420  Gel  A/D  gam. 

20430  PRINT:INPUT  *  Gain  of  A  D  (1,2,4  or  8)  (2  for  +  /-5v)*;GAIN.CODE 
20440  GAIN. CODE  =  LOG(GAIN.CODE)/LOG;2) 

20450  ’ 

20460  Get  A  D  channel. 

20470  PRINT/INPUT-  Signal  A  (Front  Beam>-USE  A/D  CH  #1*;S.CHANNEL:PRINT 

20480  E.CHANNEL  =  S. CHANNEL*  1 

20490  PRINT'SIGNAL  B  (REAR  BEAM)  ON  CH«* 

20500  PRINT: PRINT: PRINT* Press  any  key  to  continue* 

20510  AS  =  INKEYS  AS  =  "  THEN  20510 
20520  RETURN 
21000  ' 

21010  '  Instruct  the  A/D  to  lake  data 
21020  ' 

21030  WAIT  STATUS. REGISTER.COMM AND. WAIT 

21040  OUT  COMMAND. REGISTER.CSAD 

21050  WAIT  STATUS. REGISTER.WRITEWAIT.WRITRWAIT 

21060  OUT  DATA  REG ISTER.G AIN. CODE 

21070  WAIT  STATUS. REGISTER, WRITE. WAIT, WRITE. WAIT 

21080  OUT  DATA.REGISTER.START.CHANNEL 

21090  WAIT  STATUS.  REGISTER,  WRITE  WAIT, WRITE  WAIT 

21100  OUT  DATA. REGISTER. END. CHANNEL 

21110  WAIT  STATUS.  REGISTER.WRITE  WAIT, WRITE  WAIT 

21120  OUT  DATA  REGISTER.DUMMY 

21130  WAIT  STATUS.  REGISTER, WRITE  WAIT, WRITE  WAIT 

21140  OUT  DATA.REGISTER.DUMMY 

21150  WAIT  STATUS.  REGISTER.COMM  AND.  WAIT 

21 160  OUT  COMMAND  REGISTER, CCLOCK 

21170  PH#  =  INT(TICKSk/256) 

21175  PL<*  =  TICKSfc-PH«f,256 
21180  PERIODH  =  PH# 

21185  PERIODL  =  PLa 

21 190  WAIT  ST  ATUS.REGISTER, WRITE  WAIT, WRITE  WAIT 

21200  OUT  DATA.REGISTER, PERIODL 

21210  WAIT  ST  ATUS.REGISTER, WRITEWAIT.WRITEWAIT 

21220  OUT  DATA.REGISTER, PERIODH 

21230  DMACOUNT  =  (NUM.CO.NV*2)-l 

21240  DMACOUNTH  =  1NT( DMACOUNT "2 56) 

21250  DMACOUNTL  =  DMACOUNT-DMACOUNTH*256 
21260  OUT  11, DMAMODE  ’  set  DMA  mode 

21270  OUT  12,0  '  clear  byte  flip-flop 

21280  OUT  BASEREG.DMABASEL  '  set  DMA  memory  base  address 
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21290  OUT  BASEREG.DMABASEH 

21300  OUT  COUNTREG.DMACOUNTL  '  set  DMA  byte  count 
21310  OUT  COUNTREG.DMACOUNTH 

21320  OUT  PAGEREG.D.MAPAGE  '  set  DMA  memory  page 
21330  OUT  10.DMACHANNEL  '  enable  DMA  channel  mask 

21340  WAIT  STATUS. REGISTER, COMMAND. WAIT 
21350  STATUS  =  1NP(STATU’S. REGISTER) 

21360  IF  (STATUS  AND  &H80)  THEN  GOTO  22000 

21370  WAIT  STATUS. REGISTER, COMMAND. WAIT 

21380  OUT  COMMAND. REGISTER, CRAD  CDMA  +■  COMMAND 

21390  WAIT  STATUS. REGISTER, COMMAND. WAIT 

21400  STATUS  =  INP(STATUS. REGISTER) 

21410  IF  (STATUS  AND  &H80)  THEN  GOTO  22000 
21420  RETURN 
22000  ’ 

22010  'Error  handling  for  A/D 
22020  ' 

22030  PRINT'FATAU  BOARD  ERROR' 

22040  PRINT'STATUS  REGISTER  VALUE  IS  *;HEXS(STATUS);'  HEXIDECIMAL' 
22050  PRINTGOSUB  22100 

22060  PRINT-ERROR  REGISTER  VALUES  ARE:' 

22070  PRINT'  BYTE  1  -  ';HEXS(ERRORI );'  HEXIDECIMAL' 

22080  PRINT'  BYTE  2  -  ■;HEXS(ERROR2);'  HEXIDECIMAL' 

22090  PRlNT:GOTO  22240 

22100  ' 

221 10  Read  the  error  register 

22120  OUT  COMMAND. REGISTER, CSTOPTEM P  =  INP(DATA.REGISTER) 

22130  WAIT  STATUS. REGISTER.COM.MAND. WAIT 
22140  OUT  COMMAND  REGISTER.CERROR 
22150  WAIT  STATUS. REGISTER, READ. WAIT 
22160  ERRORI  =  INP(DATA.REGJSTER) 

22170  WAIT  STATUS.REGISTER.READ  WAIT 
22180  ERROR2  =  1NP(DATA. REGISTER) 

22190  STOP 
22200  ' 

22210  Illegal  Status  Register. 

22220  PRINT'FATAU  ERROR  -  ILLEGAL  STATUS  REGISTER  VALUE- 

22230  PRINT'STATUS  REGISTER  VALUE  IS  ';HEXS(STATUS);'  HEXIDECIMAL' 

22240  STOP 
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Appendix  B.  PROGRAM  OILCF.FOR 


ccc - PROGRAM  OILCF - 

C 

C 

C  VERSION  6  INPUTS  N  ZEROS  AND  N' ZERO  A  MAXIMUMS 
C  AND  PRINTS  TO  FOR  EACH  BEAM 

C  NOTE;  PRINTER  IS  REQUIRED 

C 
C 

CCC  REDUCES  DIGITAL  OUTPUT  OF  PHOTODIODE  TO  A  COEFFICIENT  OF 
CCC  SKIN  FRICTION  BY  USE  OF  OIL  LUBRICTION  THEORY. 

CCC 

CCC  USES  THE  FOLLOWING  SUBROUTINES  ; 

CCC 

CCC  SMTHPD--SMOOTHES  RAW  DIGITAL  OUTPUT  OF  PHOTODIODE 

CCC  SOILFM-IDENT1FIES  FRINGES  OF  PHOTODIODE  OUTPUT 

CCC  SLUBFT-APPLES  LUBRICATION  THEORY  TO  DATA 

CCC  C LFCLC-C A LC U LAT ES  THE  WALL  SHEAR  FOR  THE  OIL  FILM 

CCC 

CCC  THIS  PROGRAM  IS  BASED  ON  SUBROUTINE  USED  BY  DR.  WESTPHAL 
CCC  AT  NASA-AMES  AND  GIVEN  IN  NASA  TM-88216  'IMPROVED 

CCC  SKIN  FRICTION  INTERFEROMETER' 

CCC 

CCC 

CCC 

CCC  THIS  PROGRAM  READS  THE  INPUT  DATA  FILE  AS  A  SEQUENTIAL 
CCC  ACCESS  FILE  (FILE  * I).  PARAMETERS  ARE  READ  IN 
CCC  FROM  THE  SCREEN.  OUTPUT  FILES  WRITTEN  ARE. 

CCC  K12-SMOOTHED  INPUT  DATA  (INTEGER  COUNTS) 

CCC  «U3  &  *  14-BEAM  A&B  ZEROES  AND  TIME  LOCATIONS 

CCC  SET  DIMENSION,  REALS,  INTEGERS,  COMMONS,  AND  DATA 
CCC 

IMPLICIT  DOUBLE  PRECISION  (A-H) 

IMPLICIT  DOUBLE  PRECISION  (O-Z) 

INTEGER  IRA WA(5000),IRAWB(5000),NZA(6),NZB(6),ICHK(6), FORM 
CHARACTER*  18  SOURCE, OUTFIL 

COMMON/FITCOM.CFNFIT,ASLFIT,TSTFIT,ZERO(200),NZERO.RMSERR 
EXTERNAL  FITERR 

DOUBLE  PRECISION  MUNOM.NOIL,MIN'US,ERRA(6),ERRB(6),DT(6),CFM(6),ERO 

i r\T/*\  T n  a  /A\  Tnn(*\ 

.  w  ,  *  “"v"/ 

DO  203,  K  =  1,6 
ERRA(K)  =  0.D0 
203  ERRb(K)  =  0.D0 
NUL  =  0 

WRITER, •)'  INPUT  NUMBER  OF  ZEROES  TO  USE  FOR  BEAM  A' 

READ(*,*)NA 

WRITE(*,*)'  INPUT  NUMBER  OF  ZEROES  TO  USE  FOR  BEAM  B' 

READ(*,*)NB 

CCC 

CCC 

CCC  OPEN  INPUT  AND  OUTPUT  FILES 

CCC 

CCC 

CCC  INPUT  FILE  #1— RAW  DIGITAL  PHOTODIODE  OUTPUT 


Appendix  B.  PROGRAM  OILCF.FOR 


104 


WRITE(Ml) 

1 1  FORM  AT(5X,'ENT£R  FILE  #1  .DATA') 

OPEN(l .FILE  =  '  '.STATUS  =  OLD  ) 

INQUIRE!  1  .NAME  =  SO  U  RCE) 

OPEN(9,FlLE  =  LPT1  '.STATUS  =  OLD) 

WRlTE(9,*)SOURCE 

ccc 

CCC  OUTPUT  FILE  *12— SMOOTHED  INPUT  DATA  (INTEGER) 

WR1TE(*,13) 

13  FORM AT(5X, 'ENTER  FILE  *12, OUTPUT') 

OPEN(l  2.F1LE  =  *  '.STATUS  =  NEW) 

CCC  OUTPUT  FILE  *13-BEAM  A  ZEROES  AND  THEIR  LOCATIONS 
WRITE(*,14) 

14  FORM AT(5X, 'ENTER  FILE  *13,OUTPUT) 

OPEN(13,FILE  =  '  '.STATUS  = 'NEW) 

CCC  _  OUTPUT  FILE  *14— BEAM  B  ZEROES  AND  THEIR  LOCATIONS 
VVRITEC.15) 

15  FORMAT(5X,'ENTER  FILE  *14, OUTPUT') 

OPEN(14,FILE=  '  'STATUS  =  NEW) 

CCC 

CCC  READ  IN  FREQUENCY  OF  DATA  TO  COMPUTE  THE  DELTA-TIME 
CCC  BETWEEN  DATA  POINTS  (FREQ  IN  HERTZ) 

WRITE(M) 

1  FORMAT(2X.ENTER  FREQUENCY  (IN  HERTZ)  OF  DATA  TO  COMPUTE  DELTA-T) 
READC, *)FREQ 
DELT=  IDO  FREQ 

CCC  READ  IN  PARAMETER  VALUES  FOR  THE  FOLLOWING: 

CCC 

CCC  PARAMETERS  FOR  THE  SMTHPD  SUBROUTINE 
CCC  NRAWMX-N UMBER  DATA  PTS,  LENGTH  OF  IRAWA  &  IRAWB  ARRAYS 
CCC  FACTOR-SMOOTHING  FACTOR,  BTWN  1  &  0  1-FLAT  CURVE  O-NO  SMOOTH 
CCC  MODEP-TYPE  OF  SMOOTHING. BTWN  0  &  10.0-CAUSAL.1-3  PT  (ALSO  * 

CCC  OF  TIMES  3  PT  SMOOTHING  IS  USED  IF  1-10) 

WRITE(V50) 

50  FORMAT(5X, ENTER  FACTOR.MODEP  ) 

READC  ,*)FACTOR.MODEP 

CCC 

CCC 

CCC  PARAMETERS  for  the  soilfm  subroutine 

CCC  MINFA&B-MIN  ARGUMENT  FOR  DATA  REDUCTION  (MIN  DIMENSION  OF  RAW  DATA) 
CCC  MAXFA&B-MAX  *  *  (MAX  DIMENSION) 

CCC  NINTA&B-MIN  *  OF  PTS  IN  ONE  INTERVAL  USED  TO  FIND  ZERO  CROSSING 
CC  LINE.  GREATER  THAN  I 

CCC  RMSFA&B-FR ACTION  OF  RMS  TO  USE  FOR  TESTING  IF  PT  IS  IN  REGION 
CCC  NEAR  ZERO 

CCC  NZMAXA&B-AMT  OF  STORAGE  FOR  ZERO(  )  IN  CALLING  PROGRAM 
WRITEC.51) 

51  FORMAT(5X, 'ENTER  MINFA.MAXFA.NINTA.RMSFA.NZM AXA  &  -B  VALUES') 
READ(*,*)MINFA.MAXFA..NTNTA.RMSFA,NZMAXA,MINFB,MAXFB.NINTB. 

&RMSFB.NZMAXB 

CCC 

CCC 

CCC  PARAMETERS  FOR  THE  SLUBFT  SUBROUTINE 

CCC  MODEA&B-SPECIFIES  IF  INITIAL  TIME  BASE  TO  BE  FIT  OR  AN  INPUT, 

CCC  MODE  =  1  TSTFP  IS  INPUT  MODE  =  2  TSTFP  FOUND  FROM  FIT 

CCC  TSTFPA  & -B  ARE  THE  TIME  BASE  FOR  FIT 

WRITEf,52) 

52  FORM  AT  (  j  X ,  ’ENTER  MODEA, MODES, TSTFFA.TSTFPB') 

READC, *)MODEA,MODEB,TSTFPA.TSTFPB 

CCC 

CCC 

CCC  PARAMETERS  FOR  THE  FRICTION  CALCULATION 
CCC  TCOIL-OIL  TEMP  IN  DEG  CELSIUS 

CCC  TC.NOM-NOMJNAL  OIL  TEMP  FOR  USE  IN  DOW  CORNING  EQN  IN  DEG  CEL 
CCC  MUNOM-NOMINAL  OIL  VISCOSITY  FOR  USE  IN  DOW  CORNING  EQN 
CCC  NOIL-OIL  INDEX  OF  REFRACTION 

CCC  WLLSER- LASER  WAVELENGTH  IN  METERS,  0.6328X10-6  FOR  HENE 
CCC  ANGIN-LASER  BEAM  INCIDENCE  ANGLE  IN  DEG 

CCC  DXBEAM-BEAM  SEPARATION  IN  METERS  FROM  FRONT  TO  REAR  BEAM 
WLLSER  =  6328D-06 

WRITE(V)  ENTER  BEAM  SPACING  IN  mm 
READC, *)DXBEAM 
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DXBEAM  =  DXBEAM.'IOOO.DO 
MUNOM  =  .0479538D0 
NOIL=  1.40 125  DO 

WRITE(*,*)’  ENTER  ANGLE  OF  INCIDENCE  IN  DEGREES' 

READ(V)ANGIN 
TCNOM  =  25. DO 

ccc 

ccc 

CCC  FLOW  PARAMETERS 

CCC  QE-DYNAM1C  PRESSURE  OF  FLOW,  IN.  WATER 
CCC  READ  QE  FROM  SCREEN 
WRITE(’J) 

3  FORMAT(5X, ’ENTER  QE(DYNAMIC  PRESS)  IN  in.  water  ) 

READ(V)QE 

C  CONVERT  QE  FROM  IN.  WATER  TO  N/M 
QE  =  QE*248.9460 1 7 1  DO 
CCC  READ  01  TEMP  FROM  SCREEN 
WRITE(*,35) 

35  FORMAT(SX, ENTER  TC01L  (OIL TEMP  IN  DEG  C)  AND  PRESSURE 
&GRAOIIBNT  (IN  N/M 
READ(*,*)TCOIL,DPDX 

CCC  NUMBER  OF  DATA  POINTS  TO  PROCESS 
WRITE(*,4) 

4  FORM  AT(5X, 'ENTER  NUMBER  OF  DATA  POINTS  TO  PROCESS’) 
READ(V)NRAWMX 

1TYP  =  1 

READ(UO)X,XVAL 
READ(1,30)Y.YVAL 
READ(1,30)Q,QVAL 
30  FORMAT(A2,F5.2) 

CCC 

CCC 

CCC  READ  IN  DIGITAL  INPUT  DATA  SIGNAL  FROM  A/D  CONVERTER 
CCC  FOR  THE  PHOTODIODE  OUTPUT 
DO  40  1  =  l.NRAWMX 
READ(1, VEND  =4000)1  RAWA(I).IRAWB(I) 

40  CONTINUE 
CCC 

CCC  DO  CALCULATIONS 
CCC 

CCC  CALL  SUBROUTINES  TO: 

CCC  l)SMOOTH  OUTPUT  OF  PHOTODIODE 

CCC  2)1. D.  FRINGES 

CCC  3)FIT  LUBRICATION  THEORY  TO  DATA 

CCC  4)CALCULATECF  BASED  ON  TYPE  OF  OIL 

CCC 
CCC 

CCC  DO  THE  FOLLOWING  SUBROUTINES  FOR  EACH  BEAM 

CCC 

CCC 

CCC  DO  SMOOTHING  FOR  BOTH  BEAM  A&B  AND  WRITE  THE  SMOOTHED 

CCC  OUTPUT  TO  FILE  #12 

CCC 

4000  CONTINUE 
CLOSE(l) 

WRITE(*, 100) 

100  FORMAT(5X,’ENTRY’) 

OvLLSMTHPD(IRAWA,NRAWMX,  FACTOR,  MODEP) 

CALL  SMTHPD(IRAWB,NRAWMX, FACTOR, MODEP) 

WRITE(*,55) 

55  FORMAT(5X,-SMTHPD  CHECK  FOR  BEAMS  A&B-NOW  WRITE  DATA  TO  #12  ) 
DO  600  1  =  l.NRAWMX 
WRITE(  1 2,*)I  RA  WA(I),IRA  WB(I) 

600  CONTINUE 
CLOSE  (12) 

OPEN(2,FlLE  =  ’SOURCEDAT’^TATUS  =  ’NEW’) 

WRlTE(2,6674)SOURCE,'OUT’ 

REWIND  2 

READ(2, 6675)0  UTFIL 
CLOSE(2) 

OPEN(3,FILE  =  OUTFIL,STATUS  =  ’NEW') 

CCC 
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CCC  CONTINUE  CALCULATIONS  FOR  ONE  BEAM  AT  A  TIME, 

CCC  FIRST  FOR  BEAM  A. 

CCC 

54  CALLSO!LFM(IRAVVA,NRAWMX,MINFA.MAXFA,N!NTA,AVGA,RMSA,RMSFA, 
&ZERO.NZMAXA.NZERO.ITYP) 

TH1CKA  =  WLLSER*REAL(NZEROA) 

C 

C  MAX  NZEROA  =  NA 
C 

IF(NZERO.GT.NA)  THEN 
L  =  0 

DO  700  K  =  NZERO-NA+  l.NZERO 
L  =  L  +  I 

700  ZERO(L)  =  ZERO(K) 

NZERO  =  NA 

ENDIF 

\VRITE(*,56) 

56  FORM AT(5X,  SOI LFM  A  CHECK-AVG, RMS, NZERO') 

WRITE(*,*)AVGA,RMSA,  NZERO 
IF(NZERO.EQ.O)  THEN 

WRITE)*,*)'  NZEROA  =  0  FOR  THIS  CASE' 

WRITE)*,*)'  ITYP  =',ITYP 

NZEROA  =  0 

NZERO B  =  0 

AERR  =  0.D0 

BERR  =  0.D0 

GO  TO  6665 

ENDIF 

NZEROA  =  NZERO 
DO  500  I  =  l.NZERO 
ZERO(I)  =  ZERO(I)*DELT 
WRITE)  1 3,*  )l  .ZERO)  I) 

500  CONTINUE 

AZEROl  =  ZERO(l) 

CALL  SLUBFT(ZERO, NZERO, MODEA.CFNFPA.ASLFPA.TSTFPA, ERROR) 
aerr  =  error 

WRITE(*.57)CFNFPA.ASLFPA.TSTFPA 

57  FORMAT(5X,  SLUBFT  CHK,  END  OF  BEAM  A  CHK-CFNFP.ASLFPA.TSTFPA', 

&,  .F12.4.1  X.F12.4.I  X.F12.4) 

CCC  SECOND  FOR  BEAM  B-  REAR  BEAM 

CALL  S01LFM(IRAWB.NRAWMX.MINFB.MAXFB.NINTBAVGB,RMSB,RMSFB, 
&ZERO.NZMAXB, NZERO. ITYP) 

THICKB  =  WLLSER*  NZERO  B 
C 

C  MAX  NZEROB  =  NB 
C 

IF(NZERO.GT.NB)  THEN 
L  =  0 

DO  710  K  =  NZERO-NB  +  l.NZERO 
L  =  L  *  1 

710  ZERO(L)  =  ZERO(K) 

NZERO  =  NB 
ENDIF 
WRITE)*  .591 

59  FORMAT(5X,  SOILFM  B  CHECK-AVG, RMS.NZERO  ) 
WRITE)*,*)AVGB.RMSB,NZERO 
IF(NZERO.EQ.O)  THEN 
WRITE)*,*)'  NZEROB  =  0  FOR  THIS  CASE' 

WRITE)* ,*)  ITYP  =  '.ITYP 

NZEROA  =  0 

NZEROB  *  0 

BERR  =  0.D0 

GO  TO  6665 

ENDIF 

NZEROB  =  NZERO 

DO  7010  1  =  l.NZERO 

ZERO(I)  =  (ZERO(I)*DELT)  +  .05D0 

WRITE)  14,*)!  .ZERO)  1) 

7010  CONTINUE 

BZEROl  =  ZERO(I) 

CALL  SLUBFT(ZERO,NZERO, MODES, CFNFPB.ASLFPB.TSTFPB, ERROR) 

BERR  =  ERROR 


Appendix  B.  PROGRAM  OILCF.FOR 


WRITER, 60) 

60  FORMAT(5X,'SLUBFT  B) 

WRITE/*, 70 1  )CFNFPB,AS  LFPB, TSTFPB 

701  FORMAT/SX.SLUBFT  CHK,  END  OF  BEAM  B  CHK-CFNFP,ASLFPB,TSTFPB\ 
&/,F12.4,lX,Fl2.4,lX,Fl2.4) 

11=0 

AVGTFB  =  ((TSTFPB  +  ZERO(NZEROB))-(TSTFPB  +  ZERO/l  )))/2D0 
3000  11  =  11  +  1 

TFCHKB  =TSTFPB  +  ZERO(II) 

IF(TFCHKB.EQ.AVGTFB)  GOTO  3200 
IFfTFCHKB.LT. AVGTFB)  GOTO  3000 
TFBPCT  =  ( AVGTFB-/TSTFPB  +  ZERO(ll)))/AVGTFB 
ENBAVG  =  (CFNFPB-FLOAT(lI))-TFBPCT 
GOTO  3500 

3200  ENBAVG  =  CFNFPB-FLOAT(II) 

3500  WRITE(*,3550)AVGTFB.TFCHKB,TFBPCT.ENBAVG 

3550  FOR\lAT(2X,'AVG  EFFECTIVE  TIME, CHECK  FOR  AVG  EFFTIME',/. 

A'PCT  DIFFERENCE  FOR  AVG  EFF  TIME,  AVG  EFF  FRINGE  NUMBER',/, 
&F12.5,2X.F12.5,2X,F8.4,2X,F8.4) 

CCC  CALCULATE  C(X2)-C(X1)  FOR  THE  SHEAR  CALCULATION 
DELNT  =  1  DO  (ASLFPB-ASLFPA) 

WRITE/*,  1080)DELNT,aerr,berr 

1080  FORMAT/2X,  DELTA  CX  CHK-DELNT,aerr,berr'/3(2X,F12.4)) 

CCC  DO  ROUTINE  TO  CALCULATE  THE  SHEAR 

CCC  THEN  USE  SHEAR  TO  FIND  FRICTION  COEFFICIENT 

CALL  CFCLC(TCO!L,TCNOM,MUNOM,N'OIL,WLLS£R,ANG!N,DXBEAM, 
CDPDX.ENBAVG.DELNT.TAU.TAUBAR.EPSl.IER) 

CF  =  TAUBAR/QE 
CCC 
CCC 

CCC  WRITE  OUTPUT  DATA  TO  SCREEN 

CCC  OUTPUT  DATA  IS  TAUBAR,CF,QE,AVG,RMS,NZERO,CFNFP,TSTFP,ASLFP 
CCC 

4010  FORMAT(5X,TAUBAR  =  '.F16.8,/,5X,CF  =  ',F16.8,/,5X.'QE  =  \ 

&F16.8) 

WRITE(*,4010)TAUBAR,CF,QE 

WRITE(*,4020)AVGA,RMSA,NZEROA,CFNFPA,TSTFPA.ASLFPA 
WRITE(*,4030)AVGB,RMSB,NZEROB,CFNFPB, TSTFPB, AS  LFPB 
WRITE/*,*)'  HEIGHT  OF  OIL  FILM  =  '.CFNFPA* WLLSER*  1 E6,'  microns' 

6665  CONTINUE 

DENOM  =  (TSTF  PA  +  AZERO 1  +  TSTFPB  +  BZERO I  )*.5d0 

!F(!TYP.LT.6)  THEN 

CFM(ITYP)  =  CF 

ERRA(ITYP)  =  AERR 

ERRB(ITYP)  =  BERR 

TOA(ITYP)  =  TSTFPA 

TOB(ITYP)  =  TSTFPB 

NZA(ITYP)  =  N  ZERO  A 

NZB(ITYP)  =  NZEROB 

DT(ITYP)  =  ABS(TSTFPA-TSTFPB)/DENOM 

ENDIF 

IF/ITYP.EQ.6)  THEN 
ERRA/6)  =  AERR 
ERRB(6)  =  BERR 
T0A(6)  =  TSTFPA 
T0B(6)  =  TSTFPB 
NZA(6)  =  NZEROA 
NZB(6)  =  NZEROB 

DT(6)  =  ABS(TSTFPA-TSTFPB)/DENOM 
CFM(6)  =  CF 

CFAVG  =  (CFM(2)  +  CFM(3)  +  CFM(1)  +  CFM{4)+  CFM(5)+  CFM(6))  6.D0 
PLUS  =  DMAX1(CFM(2),CFM(3),CFM(1),CFM(4),CFM(5),CFM(6))  -  CFAVG 
MINUS  =  CFAVG  -  DMIN1(CFM(2),CFM(3),CFM(1),CFM(4),CFM(5),CFM(6)) 
DO  199  U  =  1,6 

199  EROOT(IJ)  =  SQRT(ERRA(U)*ERRA(IJ)  +  ERRB(U)*ERRB(IJ)  +  DT(IJ)* 
1DT(IJ)) 

WRJTE(*,6666) 

DO  201  IJ  =  1,6 

IF(1J.EQ.1)ASS1GN  6667  TO  FORM 
1F( I J.EQ.2) ASSIGN  6668  TO  FORM 
1F(IJ.EQ.3)ASSIGN  6669  TO  FORM 
IF(U.EQ.4)ASSIGN  6670  TO  FORM 


Appendix  B.  PROGRAM  OILCF.FOR 


108 


IF(IJ.£Q.5)ASSIGN  6671  TO  FORM 
IF(IJ.EQ.6)ASS1GN  6672  TO  FORM 

WRITE(3,FORM)NZA(U),NZB(IJ).CFM(U).ERRA(lJ).ERRB(IJ),DT(m 

201  WRITE(*,FORM)NZA(I  J),NZB(lJ)>CFM(lJ),ERRA(li),ERRB(IJ),DT(lJ) 
WRITE(*,6673)CFAVG,  PLUS,  MINUS 

WR!TE(9.6666) 

DO  202  IJ  =  1.6 

IF(IJ.EQ.t)ASSIGN  6667  TO  FORM 
1F(1J.EQ.2)ASS1GN  6668  TO  FORM 
IF(!J.EQ.3)ASSIGN  6669  TO  FORM 
IF(IJ.EQ.4)ASSIGN  6670  TO  FORM 
IF(1J.EQ.5)ASSIGN  6671  TO  FORM 
IF(IJ.EQ.61ASSIGN  6672  TO  FORM 

202  WRITE(9.FORM)NZA(IJ)1NZB(IJ),CFM(IJ).ERRA(U).ERRB(IJ),DT([J).TO 
1  A(IJ).T0B(1J) 

WR|TE(9,*)'  ' 

C  WRITE(9,*)'  HEIGHT  OF  OIL  FILM  =  '.CFNFPA*  WLLSER*  1 E6,'  microns' 

C 

C  ***  check  for  bad  cf  values 
C 

DO  204  IJ  =  1.6 

204  ICHK(U)  =  0 
DO  205  IJ  =  1,6 
IF(NZBdJYEQ.O)  ICHK(IJ)  =  1 
lF(dtdJ).GT.  02)  ICHK(IJ)  =  I 
1F(ERRA(IJ).GT..004)  ICHK(IJ)  =  1 
IF(ERRB(IJ).GT..004)  ICHK(IJ)  =  1 
IF(NZAdJ).LT.lO)  ICHK(IJ)  =  1 
IF(NZB(IJ).LT.10)  ICHK(IJ)  =  1 

205  IFdCHK(IJ).EQ.t)  SUL  =  NUL  *  1 
GO  TO  5000 

ENDIF 

1TYP  =  ITYP  +  1 
GO  TO  54 
CC 

CC  . WRITE  TO  OUT  FILE . 

CC 

5000  CLOSE  (3) 

CLOSE  (13) 

CLOSE  (14) 

6666  FORMAT(24X.  CF\IOX.  ERRA  .7X,  ERRB'.SX.  DT  ') 

6667  FORM AT(  1  X,12,7',i2,5X.O  = -.F12.8.3F1 1.5.2F16.5) 

6668  FORMAT) IX.  12,T,i2,5X,  M  =  \F12.8 JF1 1.5.2F16.5) 

6669  FORMAT(lX,I2,'/',i2,5X,'P  =  '.F12.8JF1 1.5.2F16.5) 

6670  FORMAT(  1  X,I2,7',i2,5X,'0NEG  =  '.F12.8 JFl  1.5.2F16.5) 

6671  FORM  AT(1  X.I2,7.i2.5X.'MNEG  =  .F12.8.3F1 1 .5.2F16.5) 

6672  FORMAT)  1  X.I2,7'.i2,5X.'PNEG  =  \F12.8,3F1 1.5.2F16.5) 

6673  FOR.MAT(20x,'CF-AVG  =  •.F12.8/20X,'.MAX-PLLS  =  +  '.F9.5 /20X.  MAX-MINU 
*S  =  -,F9.5) 

6999  FORMAT(F12.8,4Fl  1.5) 

6674  FORM ATf A18.tl3,a3) 

6675  FORMAT(8X,A10) 

4030  FORMAT(2X,AVGB  =  '.Flb.S.UX.RMSB  =  ',F16.8,:^X.'NZEROB  = 

&14./.2X,  CFNFPB  =  .F16.8,;JX,TSTFPB=  '.F16.8./.2X, 

&ASLFPB  =  .F16.8) 

4020  FORM AT(2X,'AVGA  =  \F16.8./,2X.'RMSA  =  \F16.8„  .2X.NZEROA  =  '. 

&I4,  ,2X.CFNFPA  =  ’.F16.8./JX.TSTFPA  =  '.F16.8./.2X. 

&ASLFPA  =  .F16.8) 

STOP 

END 

CCC 

CCC 

CCC 

CCC  ••••••••SUBROUTINES************** 

CCC 

CCC 

SUBROUTINE  S MTH PD( I RAWDT, NRAWMX, FACTOR. MODEP) 

CCC 

CCC  SUBROUTINE  SMTHPD-SMOOTH  OUTPUT  OF  PHOTODIODE 
CCC 

IMPLICIT  DOUBLE  PRECISION  (A-H) 

IMPLICIT  DOUBLE  PRECISION  (O-Z) 
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INTEGER  IRAWDT(NRAWMX) 

IF(MODEP.GT.O)  GOTO  250 
1F( FACTOR. LE  O.)  GOTO  999 
IFfFACTOR.GT.  I.)  GOTO  999 
FACT1M  =  1. -FACTOR 
SMOOTH  =  DBLE(IRAWDTfl)) 

DO  200  1  =  l.NRAWMX 

SMOOTH  =  SMOOTH ’FACTOR  *  DBLE(IRAWDT(!))*FACT1  M 
IRAWDT(I)  =  INTfSMOOTH  -  .5) 

2C0  CONTINUE 
GOTO  999 
250  CONTINUE 

UMAX  =  NRAWMX-1 
DO  400  J  =  l,MODEP 
OLD  =  DBLE(IRAWDTfl)) 

DO  300  1  =  UMAX 

SMOOTH  =  (OLD  -  D8LE(2*lRAWDT(I)-URAWDTfI  +  l)))/4.0 
OLD  =  DBLE(IRAW'DTfl)) 

1RAWDT(I)=  INTfSMOOTH-1- .5) 

300  CONTINUE 

SMOOTH  =  fOLD-  3  O’DBLEf.N  RAWMX))  4.0 
IRAWDT(NRAWMX)  =  INTfSMOOTH  -*•  .5) 

400  CONTINUE 
999  RETURN 
END 
CCC 
CCC 
CCC 
CCC 
CCC 

SUBROUTINE  S01LFM(IRAWDT,NRAWMX.M1NF1T..MAXFIT,.NTNTMN.AVG.RMS, 
ARMSFCT.ZERO.NZEM  AX.NZERO.ITYP) 

CCC  SUBROUTINE  SOILFM1 — FRINGE  I  D.  ALGORITHM,  FINDS  ZERO  CROSSINGS 
CCC  OF  EACH  FRINGE 

IMPLICIT  DOUBLE  PRECISION  (A-H) 

IMPLICIT  DOUBLE  PRECISION  (O-Z) 

DOUBLE  PRECISION  AVG.RMS.ZEROfNZEM  AX) 

INTEGER  IRAWDTfNRAWMX) 

IFfMINFIT  LT  D  MINFIT  =  1 

IF(M  AX  FIT. GT.N  RAWMX)  MAXFIT  =  NRAWMX 

I  F(( RMSFCT.LT.  1  DOVOR  f  R.MSFCT.GT.2.D0))  GOTO  990 

1F!(NINTMN.LT.1)  OR.(N!NTMN.GT.10))  GOTO  990 

AVG  =  0  0DO 

RMS  =  O.ODO 

DO  300  1  =  MINFIT, MAXFIT 
1IND  =  l-MINFIT-*-  I 

AVG  =  AVG  *  (DBLE(IRAWDT(I)>-AVG)  DBLEfIIND) 

300  CONTINUE 

DO  305  I  =  MINFIT. MAXFIT 

JIN D  -  I-MINFIT 1 

DEVIAT  =  DBLEIIRAWDTIDVAVG 

RMS  =  RMS  *  (DEVIAT*DEVIAT-RMS)/DBLEfIIND) 

305  CONTINUE 

RMS  =  DSQRTfRMS’l 

IFfITYP  EQ-2.OR.ITYP  EQ  5)  AVG  =  AVG  -  ,1*RMS 
IFfITYP.EQ  3.0R.ITYP.EQ  6)  AVG  =  AVG  +  .I*RMS 
NZERO  =  0 
IENFL  =  0 
NIN  =0 

DO  500  I  =  MINFIT, MAXFIT 

IF(DABS(DBLE<IRAWDT(I)>-AVG).LT  (RMS'R.MSFCT))  GOTO  410 
402  CONTINUE 

IF((IENFUEQ.l)  AND.(NIN.GE.NINTMN))  GOTO  420 
405  CONTINUE 
!ENFL  =  0 
NIN  =0 
GOTO  490 
410  CONTINUE 

IF(IENFLEQ.O)  IMIN  =  I 
IMAX  =  I 

NIN  =  IMAX-IMIN  +  I 
IENFL *  I 
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IF  (l.EQ.MAXFIT)  GOTO  402 
GOTO  490 
420  CONTINUE 

IF  (ITYP.LE.3)  THEN 

IF  {DBLE(lRAWDT(lMi\)).GE.AVG)  GOTO  405 
IF  (DBL£(lRAWDT(lMAX)).l.d.AVG)  GOTO  405 
ELSE 

IF  (DBLE(IRAWDT(IMIN)).LE.AVG)  GOTO  405 

IF  (DBLE(lRAVi  DT(1MAX)).GE.AVG)  GOTO  405 

END1F 

SX  =  0.0  DO 

SY  =  O.ODO 

SXX  =  0.0D0 

SYY  =  0.000 

DO  450  JJ  =  IMIN.IMAX 
RJ)=  DBLEiJJ-lMIN) 

SX  =  SX  -  RJJ 
SXX=SXX-  RJJ'RJJ 
SY  =SY  -  DBLEnRA\VDT(JJ'i)-AVG 
SXY  =  SXY  -  (DBLEdRAW  DT( Ji))-AVG)* RJJ 
450  CONTINUE 

D  =  DBLE(N[NT*SXXSX*SX 
A  =  fSXX*SY-SX*SXY’>  D 
B  =  (DBLE(N1N'I*SXY-SX*SY)  D 
ZDUM  =  -A  B-  DBLE(IMIN) 

W  R I  TEC  .30001  B .  I M I N ,  I M  AX  ,ZD  U  M 
3000  FORMAT!  B=  .El 5.7,’  IMIN  =  \I5,  [MAX=  ’.IS.' ZERO  =  ’.El 5.7) 

IF(B.LE.0/l  GOTO  405 

IF((ZDUM.LT.DBLE(IMIN)).OR.(ZDU.M.GT.  D8L£(IMAX)))  GOTO  405 
NZERO  =  NZERO  -  1 

lF(.NZERO.GT.NZEMAX)  STOP  SCZERO:  LACK  OF  STORAGE' 

ZERO! NZERO)  =  -A,B  -  DBLE(IMIN) 

GOTO  405 

499  CONTINUE 

500  CONTINUE 
GOTO  999 

990  STOP  SCZERO  .  BAD  ARGUMENT  VALUE' 

999  RETURN 
END 

CCC 

CCC 

CCC 

CCC 

CCC  SUBROUTINE  SLUBFT  TO  FIT  DATA  TO  LINEAR  LUBRICATION  THEORY 
CCC 

SUBROUTINE  SLUBFT(ZEROP,NZEROP,MODELF,CFNFP,ASLFP,TSTFP,ERRGR) 
IMPLICIT  DOUBLE  PRECISION  (A-H) 

IMPLICIT  DOUBLE  PRECISION  (O-Z) 

DOUBLE  PRECISION  PARAM(2),DELTA(2).ZEROP(200) 

COMMON  FITCOM.'  CFNFIT  ,ASLFIT.TSTFIT.ZERO(200),NZERO,RMSERR 
EXTERNAL  FITERR 
ABSERR  =  0  0025D0 
TOL  =  0  00002DO 
ERRMAX  =  0.05DO 
ITMAX  =  100 
NZERO  =  NZEROP 
DO  5  I  =  I. NZERO 
5  ZERO(I)  =  ZEROP(I) 

CFNFIT  =  DBLE(NZER01 
IF(MODELF.EQ.l)  TSTFIT  =  TSTFP 
IFIMODELF.EQ.2)  TSTFIT  =  O.ODO 
DELTAI1 )  =  1.0D0 
DELTA(2)  =  ZERO(2)-ZERO(I) 

PA  RAM(I)  =  CFNFIT 

PARAM(2)  =  TSTFIT 

OLDERR  =  FITERR(MODELF.PARAM) 

ITNO  =  0 

15  ITNO  =  ITNO  *  1 

C  ALL  GRIDLS(MODELF, PARAM, DELTA, DUM, FITERR) 
\VR1TEC,2000)!TNO,DUM 

2000  FORMATCITERATION  NO.  =  ’.15,  ERROR  = '.F12.5) 

CFNFIT  =  PARAMO) 
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TSTF1T  =  PARA\1(2) 

IF(DUM.LT.ABSERR.and.dabs(dum-olderr).It.lol)GOTO  20 
IF(DABS(DUM-OLDERR).LT.TOL)  GOTO  20 
IF((lTNO.GT.ITMAX)  AND.(DUM.UT.ERRMAX))  GOTO  20 
!F(ITNO.GT.lTM  AX)GOTO  990 
OLDERR  =  DUM 
GOTO  15 

20  CFNFP  =  CFNF1T 
ASLFP  =  ASLFIT 
TSTFP  =  TSTFIT 
ERROR  =  DC M 
GOTO  999 
990  CONTINUE 

STOP  SCCF1T:  NO  CONVERGENCE’ 

999  RETURN 
END 

DOUBLE  PRECISION  FUNCTION  FITERR(MODELF.PARAM) 

IMPLICIT  DOUBLE  PRECISION  (A-H) 

IMPLICIT  DOUBLE  PRECISION  (O-Z) 

DOUBLE  PRECISION  PARAMI2) 

COMMON  FITCOM  CFNFIT,ASLF!T,TSTFIT,ZERO(200),NZERO,RMSERR 
CFN  FIT  =  PARAMO) 

IF(MODELF.EQ.2)  TSTFIT  =  PARAM(2) 

SC  =  0. 

SCC  =  0. 

DO  I  I  =  l.NZERO 

Cl  =  (CFN  FIT-D  BLEU))*  (ZERO(I)  -  TSTFIT) 

SC  =  SC-(CI-SC)  DBLE(I) 

SCC  =  SCC-(C1*CI-SCC)  DBLE(I) 

1  CONTINUE 

RMSERR  =  DSQRT(SCC-SC'SC) 

ASLFIT  =  SC 

FITERR  =  RMSERR,1  ASLFIT 
RETURN 
END 
CCC 

CCC  SUBROUTINE  GRIDLS 
CCC 

SUBROUTINE  GRIDLS(NTERMS, A, DELTAA, CHISQR, FCHISQ) 

IMPLICIT  DOUBLE  PRECISION  (A-H) 

IMPLICIT  DOUBLE  PRECISION  (O-Z) 

DIMENSION  A(NTERMS),DELTAA(NTERMS) 

CHISQR  =  O.ODO 

DO  90  J  =  l.NTERMS 

CHISQ1  =  FCHISQ(NTERMS.A) 

FN  =  O.ODO 
DELTA  =  DELTAA(J) 

41  A(J)  =  A(J)-  DELTA 

CHISQ2  =  FCHISQ(NTERMS^A) 

IF(CHISQ1-CH!SQ2)  51.41,61 
51  DELTA  =  -DELTA 
A(J)=  A(J)~  DELTA 
SAVE  =  CHISQ1 
CHISQ1  =  CHISQ2 
CHISQ2  =  SAVE 
61  FN  =  FN  -  1. 

A(J)  =  A(J) *  DELTA 

CHISQ3  =  FCHISQ(NTERMS.A) 

IF(CHISQ3-CHISQ2)  71,81.81 
71  CHISQ1  =  CHISQ2 
CHISQ2  =  CHISQ3 
GOTO  61 

81  DELTA  =  DELTA*(LDO(1.DO  (CHISQI-CHISQ2),(CHISQ3-CHISQ2))4-0.5D0) 
A(J)  =  AfJ)-DELTA 

84  DELTAA(J)=  DELTAA(J)*FN/3.D0 
90  CONTINUE 

CHISQR  =  FCHISQ(NTER.MSA) 

RETURN 

END 

CCC 

CCC  SUBROUTINE  CFCLC-TO  CALCULATE  THE  FRICTION  COEFF 
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SUBROUTINE  CFCLC(TCOIL,TCNOM,MUNOM,NOIL,WLLSER,ANGIN,DXBEAM, 
&DPDX.ENBAVG,DELNT.TAU.TAUBAR,EPSI,IER) 

IMPLICIT  DOUBLE  PRECISION  (A-H) 

IMPLICIT  DOUBLE  PRECISION  (O-Z) 

DOUBLE  PRECISION  Pl.MUOlL.OPOIL.MUNOM.NOIL 
PARAMETER  (PI  =  3.U15926D0) 

INTEGER  1ER 
IER  =  1 

IF((TCOlL.LT.0.D0).OR.(TCOIL.GT.  100.D0))  GOTO  990 
IF((TCNOM.LT.-206.D0).OR.(TCNOM.GT.200.D0))  GOTO  990 
IF((MUNOM.LT.O.DO).OR.(MUNOM.GT.  1  ,D6))  GOTO  990 
lF((\VLLSER.LT.0.D0).OR.(WLLSER.GT.10.D-6))  GOTO  990 
IF((NOIL.LT  O.DO).OR.(NOI  L.GT.5.D0))  GOTO  990 
1F((ANGIN.LT.O.DO).OR.(ANGIN.GT.90.DO))  GOTO  990 
1F((DXBEAM.LT.O.DO).OR.(DXBEAM.GT.1  .D01)  GOTO  990 
IF((DELNT.LT.-LD-6).OR.(DELNT.GT.i.D6))  GOTO  990 
MUOIL  =  MUNOM.IDEXP(O.OU6DO*ITCOIL-TCNOM))) 

ANGOIL  =  DASIN(DSINIPI*  ANGIN  180.D01  NOIL) 

OPOIL  =  WLLSER  (2.DO’NOIL’DCOS(ANGOIL)) 

TAU  =  M  UOI L*  DX  BEAM’ DELNT,  OPOIL 
EPSI  =  (OPOIL*  DPDX’ENBAVG)TAU 
TAUBAR  =  TAU  (1D0-EPSI) 

GOTO  999 
990  IER  =  0 
999  RETURN 
END 
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IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

CHARACTER*!  PT,ANG,T,END*4,OUTFIL(9)*17,REDFIL<9)*17,END2*4,PRE*7, 
I  PRE2*7 

DIMENSION  CF00(6),CF45(6),CF60f6),C0(3).C45(3).C60(3),DC0(3),DC45( 
l3),DC60(3),DCX0(3),DCX45(3).DCX60(3).DCZO(3).DCZA5(3),DCZ60(3),HTO 
2(3),HT45(3),HT60(3),AXWT(9),AZ\VT(9),CM(9),COSA(9),COSG(9),SINA(9), 
3Q(9),PX(9).PZ(9),QX(9),OZ(9),TX(9)ITZ(9),AMAT(9.10),DUM(9,10) 

OPEN(l.FILE=  DATFIL.DAT  .STATUS  = 'NEW  ) 

1000  FOR.VlAT(A2) 

1010  FORM  AT(1  X.A7.3A2.A4) 

1020  FORMAT(A181 

1030  FORMAT(lX, POINT  NUMBER'.A3,':'/1  IX, 'NO  = '.!3,5X,'N45  =',13.5X,'N6 
10  =  13/1 IX. MAGNITUDE  = '.F11.7.10X, ANGLE  = '.F7.1,' degrees',  19X,' 

2-  /-■,F7.1,'%',17X,'i  /-',F7.i,'  degrees'/) 

1035  FORMATflOX.'CXBAR  =  *,F;  1.7.20X.CZBAR  =  .Fl  1.7/9X.DCXBAR  =  '.FI 

1 1 .7, '%  19X.DCZBAR  =  ,FU.7,'%713X,'SIGCX  =  '.F11.7,'%',20X.SIG 
2CZ  =  -,F1 1.7.  %  ) 

1036  FORMAT(10X.'DCX0('.II,')  =  ',FU.7,'%',17X,DCZ0('.ll.')  =  '.Fll 

1.7. '%') 

1037  FORMAT(10X.'DCX45(',I1 ,’)  =  ',F1 1.7.'%'.17X.'DCZ45(Ml.')  =  '.Fll 

1.7, '%') 

1038  FORMAT(10X,'DCX60('.U ,’)  =  ,F1 1.7.'%'.17X.DCZ60(',I1 .')  =  '.Fll 

1.7. '%') 

1039  FORMAT(') 

1040  FORMAT(18X,Fl  1.7JF1 1.5) 

1046  FORMAT(  1 3X,  C0(*,1 1 ,')  =  '.Fl  1.7,17X.  DC0('.!1,')  =',F11.7) 

1047  FORMAT(13X,  C45(',1 1 ,')  =  ,F)  1.7.17X,  DC45('.!1.  )  =  '.Fl  1.7) 

1048  FORMAT(13X,'C60(  .1 1 . )  =  '.Fl  1 .7.1 7X.'DC60(',1 1 ,')  =  \F11.7) 

1050  FORMATflX,  DIFF  =  '.Fl 2.8,  %  ) 

1 1 1 1  FORMAT(5X,'l  ',7X.'2'.7X.'3',7X.'4',7X.'5',7X,'6',7X,'7' 

1 ,7X,  8',7X.'9',7X.'10  ) 

1112  FOR.MAT(10D8.1) 

C 

C  INITIALIZE  VARIABLES  ••••• 

C 

R2D  =  180.  (4.*ATAN(1 .)) 

PI  =  4.*ATAN(1.) 

NO  =  0 
N45  =  0 
N60  =  0 
S1GCX  =  0 
SIGCZ  =  0 
DELX0  =  0 
DELZ0  =  0 
DELX45  =  0 
DELZ45  =  0 
DELX60  =  0 
DELZ60  =  0 
DO  301  K  =  U 
C0(K)  =  0 
DC0(K)  =  0 
C45(K)  =  0 
DC45(K)  =  0 
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C60(K)  =  0 
301  DC60(K)  =  0 

COS4S  =  COS(PI*.25) 

S1N45  =  COS45 
COS60  =  COS(PI/3.) 

SIN60  =  S1N(PI  3.) 

BETA  =  3.*  PI, '180. 

*MM  GENERATE  NINE  FILENAMES . 

WRITE!*, *)’  ENTER  POINT  NUMBER’ 

READ!*, 10001  PT 

WRITE(*,*)'  ENTER  INITIAL  GUESS  FOR  ALPHA’ 

READ(*,*)ALPHA 
ALPHA  =  ALPHA/R2D 
END  =  ’.OUT’ 

END2  =  ’.RED’ 

PRE=  ’C:  OUT>' 

PRE2  =  C:,RED\' 

K  =  1 

ANG  =  00’ 

1  DO  10  I  =  1,3 
!F(I.EQ.1)T  =  Tl’ 

IFfl.EQ  2)T  =  T2' 

IF(I.EQ.3)T  =  T3’ 

WRITE!  1 ,101 0)PRE2,PT,ANG,T,END2 
10  WRITE(1 ,101 0)PRE.PT,ANG,T,END 
IF(K  EQ  1)THEN 
ANG  =  45 
ELSE 

ANG  =  '60' 

ENDIF 
K  =  K  1 

IF(K.LT  4)  GO  TO  1 
REWIND  1 
DO  20  I  =  1,9 
READ(t,1020)  REDFIL(I) 

20  READ)  1 ,1020)  OUTFIL(I) 

CLOSE  (1) 

. READ  IN  CF  VALUES  ••••• 

DO  40  I  =  U 
10  =  0 

OPENfl  .FILE  =  REDFIL(I)STATUS  =  'OLD’.IOSTAT  =  ION) 

IF(ION.NE.O)  THEN 

WRITE!*,*)  REDFILE’.I,’  NOT  FOUND’ 

ELSE 

READ!  1  ,*)HT0(I) 

ENDIF 

CLOSE(l) 

J  =  1 

OPENfl, FILE  =  OUTFlL!I)STATUS  =  ’OLDMOST  AT  =  IONUM) 
IF(ION’UM.NE-O)  GO  TO  39 

29  10  =  10  ♦  I 

30  READ(l,lO40)CF00(l0).ERRA,ERRB,DT 
J  =  J  -  1 

I  FfERRA.GT.. 004. OR.ERRB.GT..004.OR.ERRA.LT.. 0000  LOR.ERRB.LT.. 0000 
1 1.0R.DT.GT..00001)  THEN 
IF(J.LT.7)  THEN 
GOTO  30 
ENDIF 
ENDIF 

IF(J.LT.7)  GOTO  29 
IF(I0.LT.2)  GOTO  39 
NO  =  NO  *  1 
HT0(N0)  =  HT0(I) 

DO  35  JJ  =  1 .10 

CF00(JJ)  =  CF00(J J)*(  1 .  -  COS(-ALPHA)*HT0(N0)) 

35  C0(N0)  =  C0(N0)  +  CF00(JJ) 

C0(N0)  =  C0(N0),  10 

FACTOR  =  (l.'HT0(N0)*COS(- ALPHA)) 
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DO  36  JJ  =  1,10 

36  DCO(NO)  =  DCO(NO)  +  (CFOO(JJ'  -  C0(N0))**2 
DCO(NO)  =  SQRT(DCO(NOV(IO-1)) 

39  CLOSE(l) 

40  CONTINUE 
DO  60  1  =  4,6 
145  =  0 

OP  EN(1  .FILE  =  R£DF1L(1),STATUS  =  OLD.JOSTAT  =  ION) 

IF(ION.NEO)  THEN 

VVRITE(*,*)'  REDFILE',1,'  NOT  FOUND' 

ELSE 

READ(1  ,*)HT45(l-3) 

ENDIF 

CLOSE(l) 

J  =  1 

OPEN(l  .FILE  =  OUTFIL(I), STATUS  =  'OLD',10STAT  =  IONUM) 
IF(IONU.M.NE.O)  GO  TO  59 

49  145  =  145  *  1 

50  READ(1.1040)CF45(I45).ERRA,ERRB,DT 
J  =  J  +  1 

if(erra.gt..oo4.or.errb.gt,.oo4.or.erra.lt..ooooi.or.errb.lt..oooo 

1 1 .0 R.DT.GT..0000 1 )  TH EN 
IF(J.LT.7)  THEN 
GOTO  50 
ENDIF 
ENDIF 

IF(J.LT.7)  GOTO  49 
IF(I45.LT.2)  GOTO  59 
N45  =  N45  +  1 
HT45(N45)  =  HT45(I-3) 

DO  55  JJ  =  1,145 

CF45(JJ)  =  CF45(JJ)*(1.  -  COS(Pl*.25-ALPHA)*HT45(N45)) 

55  C45(N45)  =  C45(N45)  +  CF45(JJ) 

C45(N45)  =  C45(N45),  145 

DO  56  JJ  =  1,145 

56  DC45(N45)  =  DC45(N45)  *  (CF45(JJ)  -  C45(N45))**2 
DC45(N45)  =  SQRT(DC45(N45)  (145-1)) 

59  CLOSE(  1) 

60  CONTINUE 
DO  80  I  =  7,9 
160  =  0 

OPEN(1,FILE=REDFIL(I)STATUS  =  'OLD',IOSTAT=  ION) 

IF(ION.NE.O)  THEN 

WR1TE(* ,*)■  REDFILE',1.'  NOT  FOUND' 

ELSE 

READ(1  ,*)HT60(I-6) 

ENDIF 

CLOSE(l) 

J  =  1 

OPEN(1  .FILE  =  OUTFlL(I)STATUS  =  OLD'.IOSTAT  =  IONUM) 
IF(IONUM.NE-O)  GO  TO  79 

69  160  =  160  +  1 

70  READ(  1 ,1040 )CF60(I60),ERRA,ERRB,DT 
J  =  J  +  1 

I  F(ERRA.GT.. 004.0  R.ERRB.GT.. 004.OR.ERRA.LT.. 0000  LOR.ERRB.LT..OOOO 
1 1. OR. DT.GT.. 00001)  THEN 
IF(J.LT.7)  THEN 
GOTO  70 
ENDIF 
ENDIF 

!F(J.LT.7)  GOTO  69 
IF(I60.LT.2)  THEN 
160  =  160-1 
GOTO  79 
ENDIF 

N60  =  N60  +  1 
HT60(N60)  =  HT60(I-6) 

DO  75  JJ  =  1,160 

CF60(JJ)  =  CF60(JJ)*(1.  -  COS(PI/3.-ALPHA)*HT60(N60)) 

75  C60(N60)  =  C60(N60)  +  CF60(JJ) 

C60(N60)  =  C60(N60),REAL!160) 

DO  76  JJ  =  1,160 
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76  DC60(N60)  =  DC60(N60)  +  (CF60(JJ)  -  C60(N60))**2 
DC60(N60)  =  SQRT(DC60(N60), (160-1)) 

79  CLOSE(l ) 

80  CONTINUE 

COMPUTE  MAGNITUDE  AND  DIRECTION 

OLDA  =  ALPHA 
NTOT  =  NO  +  N45  +  N60 

150  CONTINUE 

CALL  LSQ(NO.N45,N£3.CO,C45,C63,CX.CZ,AEL.BEL.DEL.DENOM) 

CXBAR  =  CX 
CZBAR  =  CZ 
ALPHA  =  ATAN(CZ/CX) 

WRITEC, *)'  ALPHA  =  ',ALPHA*R2D 
IF(ABS(ALPHA-OLDA).GT..05,R2D)  THEN 
DO  147  JJ  =  1.N0 

CO(JJ)  =  CO(JJ)*(1.-HTO(JJ)*COS(-ALPHA))/(1.-HTO(JJ)*COS(OLDA)) 

147  FACTOR  =  (l.-HT0(JJ)*COS(-ALPHA)) 

147  WRITEC .•)'  CO(JJ), FACTOR  , CO(JJ), FACTOR 
DO  148  JJ  =  1.N45 

C45(JJ)  =  C45(JJ)*(l.-HT45(JJ)*COS(Pl*.25-ALPHA))/(l.-HT45(JJ)*COS 
l(PI*.25-OLDA)) 

148  FACTOR  =  l.-HT45IJJ)*COS(Pl*.25-ALPHA) 

148  WRITEC, •)’  C45(JJ).FACTOR',C45(JJ), FACTOR 
DO  149  JJ  =  1.N60 

C60(JJ)  =  C60(JJ)*(l.-HT60(JJ)*COS(PI/3.-ALPHA));(l.-HT60(JJ)*COS( 

1  PI '3. -OLDA)) 

149  FACTOR  =  (l.-HT60(JJ)*COS(PI  3.-ALPHA)) 

149  WRITEC ,*)■  C60(JJ),FACTOR  ,C60(JJ), FACTOR 
OLDA  =  ALPHA 

GO  TO  150 
ENDIF 

DO  151  JJ  =  I, NO 

151  DCO(JJ)  =  SQRT(DC0(JJ)*DC0(JJ)  *  ((l.-l./(l.-TAN(-ALPHA)*TAN(BETA)) 
1)*C0(JJ))**2) 

DO  160  JJ  =  I, NO 

SIGCZ  =  SIGCZ  +  (C0(JJ)*TAN(ALPHA)-CZBAR)*,2 
160SIGCX  =  SIGCX  +  (C0(JJ)  -  CXBAR)**2 
DO  161  JJ  =  1.N45 

161  DC45(JJ)  =  SQRT(DC45(JJ)*DC45(JJ)  +  ((l.-l./(I.-TAN(PC.25-ALPHA)*T 
1  AN(BETA)))*C45(JJ))**2) 

DO  170  JJ  =  1.N45 

SIGCX  =  SIGCX  +  (C45(JJ)*COS(ALPHA)  COS(PIV25-ALPHA)  -  CXBAR)»*2 

170  SIGCZ  =  SIGCZ  +  (C45(JJ)*SIN(ALPHA)/COS(PI*.25-ALPHA)  -  CZBAR)**2 
DO  171  JJ  =  1.N60 

171  DC60IJJ)  =  SQRT(DC60(JJ)*DC60(JJ)  +  ((l.-l.,(l.-TAN(PI/3.-ALPHA)*TA 
1  N(BETA)))*C60(JJ))**2) 

DO  180  JJ  =  1.N60 

SIGCX  =  SIGCX  +  (C60(JJ)*COS(ALPHA)  COS(PI/3. -ALPHA)  -  CXBAR)**2 
180  SIGCZ  =  SIGCZ  +  (C60(JJ)*SIN(ALPHA)/COS(PI/3. -ALPHA)  -  CZBAR)**2 
SIGCX  =  SIGCX/ NTOT 
SIGCZ  =  SIGCZ/NTOT 

.....  ITERATE  TO  GET  UNCERTAINTIES  ••••• 

DO  190  K  =  1.N0 
C0(K)  =  C0(K)  +  DC0(K) 

CALL  LSQ(N0,N45,N60,C0,C45,C60,CX,CZvAEL,BEL,DEL,DENOM) 

C0(K)  =  C0(K)  -  DC0(K) 

DCX0(K)  =  (CXBAR  -  CX)"2 
190  DCZ0(K)  =  (CZBAR  -  CZ)**2 
DO  200  K  =  1.N45 
C45(K)  =  C45(K)  +  DC45(K) 

CALL  LSQ(N0,N45,N60.C0,C45.C60,CX,CZ,AEL,BEL,DEL,DENOM) 

C45(K)  =  C45(K)  -  DC45(K) 

DCX45(K)  =  (CXBAR  -  CX)**2 
C  WRITEC,*)' CX.CXBAR  = '.CX, CXBAR 
200  DCZ45(K)  =  (CZBAR  -  CZ)**2 
DO  210  K  =  1.N60 
C60(K)  =  C60(K)  +  DC60(K) 

CALL  LSQ(N0,N45,N60,C0,C45,C60,CX,CZ^EL,BEL,DEL.DENOM) 
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C60(K)  =  C60(K)  -  DC60(K) 

DCX60(K)  =  (CXBAR  -  CX)**2 
210  DCZ60(K)  =  (CZBAR  -  CZ)**2 
C 

c  •**••  COMPUTE  FINAL  UNCERTAINTY  ON  CX  AND  CZ  ••••• 

C 

DO  720  K  =  I  ,N0 

DELXO  =  DELXO  +  DCXO(K) 

220  DELZO  =  DELZO  +  DCZO(K) 

DO  230  K  =  !  ,N45 

DELX45  =  DELX45  +  DCX45(K) 

230  DELZ45  =  DELZ45  +  DCZ45(K) 

DO  240  K  =  1.N60 

DELX60  =  DELX60  +  DCX60(K) 

240  DELZ60  =  DELZ60  +  DCZ60(K) 

C 

DCXBAR  =  SQRT(SIGCX  +  DELXO  DELX45  +  DELX60) 

DCZBAR  =  SQRT(SIGCZ  +  DELZO  +  DELZ45  -  DELZ60) 

C 

C  •••••  COMPUTE  MAGNITUDE  AND  DIRECTION  AND  PRINT  OUTPUT  ••••• 

C 

CFMAG  =  SQRT(CXBAR*CXBAR  -  CZBAR'CZBAR) 

ANGLE  =  ATAN(CZBAR  CXBAR)*180./P1 

DADX  =  CZBAR  (CZBAR'CZBAR  *  CXBAR'CXBAR) 

DADZ  =  CXBAR  (CXBAR'CXBAR  -  CZBAR'CZBAR) 

DCANG  =  SQRT( DADX* DADX* DCXBAR* DCXBAR  -t-  DADZ* DADZ* DCZBAR* DCZBAR) 
DCANG  =  DCANG*  180.  PI 

DCMDCX  =  CXBAR’CXBAR  (CFMAG’CFMAG) 

DCMDCZ  =  CZBAR'CZBAR  (CFMAG'CFM AG) 

DCM  =  SQRT(DCMDCX*  DCXBAR*  DCXBAR  +  DCMDCZ*  DCZBAR*  DCZBAR) 

DCM  =  DCM /CFMAG*  100. 

DCXBAR  =  DCXBAR  CXBAR'100. 

DCZBAR  =  DCZBAR  CZBAR*  100. 

SIGCX  =  SQRT(SIGCX)  CXBAR'100. 

SIGCZ  =  SQRT(SIGCZ)<CZBAR*  100. 

OPEN(2,FILE=  LPTU, STATUS  =  OLD) 

WRIT  E(*,1030)PT.N0,N45,N60,CFMAG,  ANGLE,  DCM, DCANG 
WR1TE(*.1035)CXBAR,CZBAR.DCX8AR,DCZBAR,SIGCX,SIGCZ 
C  DO  250  JJ  =  U 

C  DCXO(JJ)  =  SQRT(DCX0(JJ))  CXBAR*100. 

C  DCZO(JJ)  =  SQRT(DCZOfJJ))  CZBAR*  100. 

C  250  WRITE(*,1036)JJ.DCX0(JJ).JJ,DCZ0(JJ) 

C  DO  260  JJ  =  1 .3 

C  DCX45(JJ)  =  SQRT(DCX45(JJ))  CXBAR'100. 

C  DCZ45(JJ)  =  SORT(DCZ45(JJ))  CZBAR*  100. 

C  260  WRITE(*,1037)JJ,DCX45(JJ),JJ,DCZ45(JJ) 

C  DO  270  JJ  =  1.3 

C  DCX60(JJ)  =  SQRT(DCX60(JJ)VCXBAR*  100. 

C  DCZ60(JJ)  =  SQRT(DCZ60(JJ))  CZBAR*  100. 

C  270  WRITE(*,1038)JJ.DCX60(JJ),JJ,DCZ60(JJ) 

WRITE(2.1030)PT.N0,N45,N60, CFMAG, ANGLE, DCM, DCANG 
C  WRITE(2.1035)CXBAR.CZBAR,DCXBAR,DCZBAR^IGCX3IGCZ 
C  DO  280  JJ  =  1,3 

C  280  \VRITE(2,lO36)JJ,DCX0(JJ),JJ,DCZO(JJ) 

C  DO  290  JJ  =  1,3 

C  290  WRITE(2.1037)JJ,DCX45(JJ),JJ,DCZ45(JJ) 

C  DO  300  JJ  =  U 

C  300  WRITE(2,1038)JJ.DCX60(JJ)JJ,DCZ60(JJ) 

C  WRITE(2,1039) 

C  DO  310  I  =  14 
C  310  WRITE(2, 1046)1, C0(!),I,DC0(1) 

C  DO  320  I  =  1,3 
C  320  WR1TE<2, 1047)1, C45(I),I,DC45(I) 

C  DO  330  I  =  14 
C  330  WRITER, 1048)I,C60(I),l,DC60(I) 

CLOSE(2) 

END 

C 

C 

SUBROUTINE  LSQ  (N0,N45,N60,CF00,CF45.CF60,CFX,CFZA,B,D,DENOM) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DIMENSION  CF00(3),CF45(3),CF60(3) 
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PI  =  4.*ATAN(1.) 

COS4S  =  COS(PI*.25) 

COS60  =  C0S(PI/3.) 

S1N60  =  SlN(PI/3.) 

A  =  NO  +  N45*.5  +  N60V25 
B  =  N45V5  +  N60*SQRT(3.)*-25 
D  =  N45V5  +  N60V75 
E  =  0. 

F  =  0. 

DO  2000  I  =  1  ,N0 
2000  E  =  E  +  CF00(i) 

DO  2010  I  =  1.N45 
E  =  E  +  CF4;(I)«COS4S 
2010  F  =  F  +  CF45(!)*COS45 
DO  2020  I  =  1.N60 
E  =  E  +  CF60(1)*COS60 
2020  F  =  F  +  CF60(I)*SIN60 
DENOM  =  A*D  -  B*B 
CFX  =  (D*E  -  B*F)/DENOM 
CFZ  =  (A*F  -  B*E)/DENOM 
RETURN 
END 
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